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Abstract 

This  research  models  and  analyzes  the  thermochemical  damage  produced  in 
Bacillus  spores  by  short,  high-temperature  exposures  as  well  the  repair  process  within 
damaged  Bacillus  spores.  Thermochemical  damage  in  spores  is  significantly  due  to 
reaction  with  water,  hydrolysis  reactions.  Applying  heat  to  the  spore  causes  absorbed  and 
chemically  bound  water  molecules  become  mobile  within  the  spore.  These  mobile  water 
molecules  react  by  hydrolysis  reactions  to  degrade  DNA  and  enzyme  molecules  in  the 
spore.  In  order  to  survive  the  thermal  inactivation,  the  spore  must  repair  the  damaged 
DNA  during  spore  germination.  The  DNA  repair  process,  as  well  as  other  germination 
functions,  is  dependent  on  reactions  catalyzed  by  enzymes  that  are  viable  after  the 
thermal  exposure.  Increased  damage  to  the  enzymes  during  thermal  inactivation,  affects 
the  rate  at  which  the  spore  DNA  is  repaired.  If  the  enzymes  are  damaged  to  such  an 
extent  that  the  DNA  repair  is  not  completed,  the  spore  is  unable  to  germinate,  or  produce 
outgrowth.  The  DNA  repair  process,  repair  enzymes,  and  outgrowth  time  influences  the 
spore’s  chance  of  survival.  Using  this  information,  a  probability  of  survival  model  was 
created  based  on  water  mobility,  hydrolysis  reactions,  an  initial  state  of  DNA,  viable 
repair  enzymes,  and  an  outgrowth  time. 
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MODELING  OF  BACILLUS  SPORES:  INACTIVATION  AND  OUTGROWTH 


I.  Introduction 


Background 

Anthrax  is  a  bacterial  infection  that  has  a  high  fatality  rate  of  89%  in  the  United 
States  that  occurs  when  Bacillus  anthracis  (B.a)  endospores  are  inhaled  into  lung  alveoli 
and  germinate  [4,  18].  Anthrax  spores  can  be  easily  prepared  from  B.a  obtained  from 
infected  animals  or  soil  [27].  The  spores  can  be  aerosolized  and  used  with  intent  to  harm 
by  causing  anthrax  disease  in  humans.  The  early  diagnosis  of  inhalational  anthrax  is  very 
difficult,  if  there  no  knowledge  of  the  anthrax  exposure.  Antibiotic  treatment  is  effective 
only  during  the  first  stage  of  this  infection,  which  is  characterized  by  nonspecific  flu-like 
symptoms.  This  stage  normally  lasts  for  several  days.  This  is  followed  by  the  final  stage, 
death  sometimes  occurs  after  only  a  few  hours  [18].  These  are  reasons  why  B.a  is  listed 
as  a  biological  warfare  agent. 

Several  incidents  have  occurred  involving  the  agent  [26].  In  1979,  anthrax  spores 
were  accidently  released  in  Sverdlovsk,  Russia.  The  release  caused  an  epidemic  of 
inhalational  anthrax  in  Russia  resulting  in  66  deaths  and  18  cases  of  exposure  to  anthrax 
in  the  United  States.  In  1993,  Aum  Shinrikyo,  a  Japanese  cult,  aerosolized  a  liquid 
suspension  of  B.a  from  the  roof  of  an  eight  story  building.  The  liquid  apparently  had  a 
low  concentration  of  spores  and  was  ineffective  since  the  B.a  strain  was  used  to  vaccinate 
animals  in  Japan.  In  2001,  anthrax  spores  were  mailed  through  the  United  States  Postal 
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Service  to  members  of  the  U.S.  Senate  and  to  national  news  agencies.  These  attacks 
resulted  in  22  combined  inhalational  and  cutaneous  cases  of  anthrax  infection.  Of  the  22 
cases,  five  Americans  died  of  inhalational  anthrax.  These  attacks  pointed  out  the 
possibilities  and  vulnerabilities  of  bioterrorism  using  B.a.  Furthermore,  these  attacks 
illustrated  the  need  for  detection,  recognition,  environmental  surveillance, 
decontamination,  and  destruction  of  anthrax. 

Previous  Research 

Much  previous  work  has  been  done  studying  and  modeling  B.a.  Researchers  at 
the  Air  Force  Institute  of  Technology  (AFIT)  have  been  researching  and  developing 
methods  to  characterize  the  effects  of  heating  B.a  spores  at  high  temperatures  for  short 
durations.  A  previous  researcher.  Captain  Emily  Knight  [19],  has  developed  a  thermal 
model  to  predict  the  effects  of  the  thermal  damage  to  the  spores.  The  thermal  model  was 
designed  to  replicate  the  laboratory  experiments  which  heat  spores  at  nearly  constant 
temperature  for  short  times.  The  thermal  model  evaluates  the  rate  of  thermal  diffusion 
throughout  the  B.a  spores  as  well  as  characterizes  the  influence  that  heating  has  on 
distributing  water  within  the  spore.  A  water  mobility  model  based  on  the  production- 
diffusion  of  water  throughout  the  spore  was  developed  to  achieve  this  goal.  Using  the 
water  mobility  model.  Captain  Knight  presented  a  probability  of  kill  model  that  relates 
repair  of  the  degradation  of  initial  DNA  information  content  to  the  decimation  of  the 
spore  population’s  protein  ‘fitness’.  This  model  was  based  on  an  assumed  Gaussian 
representation  of  protein  fitness  distribution. 
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Purpose 

The  purpose  for  this  research  is  to  develop  a  probability  of  survival  model  for  B.a 
spores  exposed  to  high  temperatures  at  short  exposure  times.  This  research  will  estimate 
the  concentration  of  water  within  the  B.a  spore  and  the  temperatures  at  which  this  water 
will  become  available  for  hydrolysis  reactions  to  occur.  Hydrolysis  is  the  reaction  of 
molecules  with  water  molecules.  Water  can  react  with  biopolymers  which  are  molecules 
such  as  DNA  and  proteins  resulting  in  deploymerization  (breaking  a  large  strand  of 
polymers  into  two  smaller  ones)  or  denaturation  (when  proteins  or  DNA  loses  their 
unique  three-dimensional  structures).  This  research  will  use  as  an  input  to  the  water 
mobility  model  mathematical  fits  of  measured  data  of  water  evolved  from  B.a  spores 
obtained  from  thermogravimetric  analysis  conducted  by  Dr.  Daniel  Felker. 

Hydrolysis  reactions  can  damage  the  spore’s  DNA  and  proteins  to  such  an  extent 
that  the  DNA  cannot  be  repaired.  The  extent  to  which  the  DNA  and  proteins  have  been 
damaged  beyond  repair  is  the  critical  threshold.  The  threshold  for  the  enzyme  damage 
influences  the  threshold  for  the  DNA.  This  is  because  the  DNA  repair  process  is 
mediated  by  the  concentration  of  viable  enzymes  (proteins  that  catalyze  specific  chemical 
reactions).  Therefore,  the  probability  of  kill  model  developed  by  Captain  Knight  [19]  will 
be  reviewed  as  well  as  different  representations  for  the  distribution  of  enzyme  fitness. 

The  spore  exist  in  a  dormant  state  where  metabolic  activity  has  stopped  and  the 
DNA  is  immobilized  in  the  spore’s  core  contributing  to  increase  the  spore’s  resistance  to 
damage.  The  spore  can  stay  in  this  dormant  state  for  long  periods.  Once  placed  in  a 
moist  and  nutrient-rich  environment,  the  spore  will  break  its  dormancy  and  begin 
germination.  During  germination,  the  spore  will  repair  damaged  DNA.  The  rate  of  this 


3 


repair  process  depends  on  the  concentration  of  viable  enzymes.  Quantitative  models  by 
Sokhansanj  et  al  [36]  exist  for  the  base  excision  repair  process.  These  models  will  be 
mimicked  to  obtain  a  qualitative  look  at  the  spore’s  DNA  repair  process.  Since  the 
enzymes  are  also  damaged  during  thermal  exposure,  the  DNA  repair  time  will  be  affected 
by  their  damage.  If  there  is  little  damage  to  the  enzymes,  then  the  repair  time  will  be 
quicker.  Whereas,  if  there  is  substantial  damage  to  the  enzymes,  the  repair  time  would 
take  longer  or  if  the  enzyme  damage  is  irreversible  the  spore  cannot  repair  the  damaged 
DNA.  For  example,  elimination  of  multiple  bases  (adenine,  guanine,  cytosine,  thymine) 
on  a  single  strand  of  the  DNA  double  helix  can  be  readily,  faithfully  repaired,  but 
completely  cleaving  the  DNA  double  helix  in  multiple  locations  may  not  be  faithfully 
repaired.  Thus,  a  spore  suffering  significant  unrepaired  DNA  damage  will  not  be  able  to 
produce  outgrowth  and,  unable  to  replicate  itself,  it  is  permanently  inactive,  i.e.  dead  for 
the  purpose  of  producing  disease. 

Many  parameters  influence  the  repair  process.  These  influences  include  the  initial 
state  of  the  spore  DNA  and  the  concentration  of  viable  enzymes.  The  concentration  of 
viable  enzymes  is  dependent  on  the  exposure  temperature  and  exposure  time. 

This  research  will  examine  the  different  influences  of  these  parameters  on  the 
repair  process  and  model  their  effect  on  the  spore’s  chances  of  survival.  Previously,  the 
probability  of  kill  model  was  based  on  an  arbitrary  estimate  of  the  enzymes’  critical 
threshold.  Unfortunately,  information  about  this  critical  threshold  is  unknown  for 
Bacillus  spores.  Therefore,  this  research  will  develop  a  probability  of  kill  model  based  on 
thermal  exposure  as  well  as  different  combinations  of  enzymes  available  for  DNA  repair 
process  and  several  outgrowth  times.  Let  outgrowth  time  be  defined  as  the  amount  of 
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time  selected  to  see  if  the  spore  produces  a  viable,  replicating  cell  during  germination. 
Therefore,  shorter  outgrowth  measurement  time  will  lead  to  a  higher  measured 
probability  of  kill  and  a  longer  measured  outgrowth  time  will  lead  to  a  shorter  measured 
probability  of  kill.  The  actual  probability  of  kill  is  estimated  by  very  long  measurement 
times.  However,  there  are  practicable  limits  to  extending  these  measurement  times. 
Producing  a  model  of  kill  probability  that  includes  damage  repair  necessary  for 
outgrowth  will  further  enhance  experimentalists’  ability  to  interpret  their  data. 

Modeling  Objectives 

•  Analyze  thermogravimetric  analysis  water  evolution  data  to  estimate  available 
water  within  the  spore  to  produce  hydrolysis  reactions  and  the  temperatures  at 
which  this  water  becomes  mobile 

•  Review  the  probability  of  kill  model  based  on  the  enzyme  threshold  killing 
mechanisms  of  B.a  spores  and  extend  to  an  enzyme  fitness  model  using  a  beta 
distribution  rather  than  Gaussian  representation 

•  Numerically  model  the  DNA  repair  mechanism  within  the  spore  and  develop  a 
probability  of  kill  model  based  on  the  DNA  critical  threshold,  DNA  repair 
process,  exposure  time,  temperature,  viable  repair  enzymes,  and  outgrowth  times. 

Overview 

The  model  development  will  be  presented  in  the  following  manner.  Chapter  II 
will  review  Captain  Knight’s  water  mobility  model.  Chapter  III  describes  the  two  types 
of  water  released  within  the  spore  due  to  heating.  The  thermogravimetric  analysis 
conducted  will  be  used  in  this  chapter  to  determine  the  temperatures  at  which  these 
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waters  are  available  within  the  spore  for  hydrolysis  reactions.  Chapter  IV  reviews  the 
probability  of  kill  model  and  extends  that  ideal  to  employ  the  more  physically  realistic 
beta  distribution.  Chapter  V  examines  the  model  of  the  spore’s  probability  of  survival 
after  the  thermal  inactivation  treatment.  Since  the  spore’s  DNA  is  damaged,  the  damage 
DNA  must  be  repaired  in  order  to  synthesize  the  enzymes  and  produce  outgrowth. 
Chapter  V  will  also  describe  and  present  results  from  the  probability  of  kill  model  based 
on  the  concentration  of  repair  enzymes  available  for  the  DNA  repair  process.  Chapter  VI 
will  give  the  conclusions  and  future  work. 
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II.  Spore  Structure  and  Dynamics 


Introduction 

The  spore’s  strong  resistance  to  heat  and  chemicals  is  primarily  due  to  the 
properties  of  multiple  layers  that  make  up  the  spore.  These  layers  protect  the  spore’s 
DNA  from  damage  by  acting  as  a  diffusion  barrier  and  by  dehydrating  the  spore’s  core 
where  the  DNA  resides.  During  the  thermal  inactivation  process,  water  within  the  spore 
becomes  mobile  and  can  react  with  molecules.  In  addition  to  pyrolysis  reactions  in  which 
bonds  that  hold  the  molecule  together  break  by  heating,  the  mobile  water  can  also  cause 
hydrolysis  damage  to  the  spore’s  DNA  and  enzymes.  Further,  the  release  and  movement 
of  molecular  water  is  important  for  these  hydrolysis  reactions. 

Spore  Structure 

The  Bacillus  spore  is  constructed  of  several  layers  including  the  exosporium, 
spore  coat,  germ  cell  wall,  outer  and  inner  membrane,  cortex,  and  core.  The  structure  and 
chemical  composition  of  the  spore  play  major  roles  in  spore  resistance.  Figure  2.1  shows 
an  illustration  of  the  major  layers  of  a  Bacillus  spore  (note  that  layers  are  not  drawn  to 
scale).  The  outermost  layer  is  the  exosporium.  The  exosporium  is  not  present  in  spores 
of  all  Bacillus  species.  It  contains  mostly  water,  but  also  carbohydrates  and  proteins  [33]. 
However,  the  exosporium  is  not  believed  to  provide  extra  protection  or  heat  resistance  for 
the  spore  [16]. 
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Beneath  the  exosporium  is  the  spore  coat  made  up  largely  of  proteins.  The  spore 
coat  protects  the  spore  from  destruction  by  lytic  enzymes  and  many  toxic  chemicals.  [33]. 
However,  heat  resistance  is  not  provided  by  the  spore  coat  [6]. 

Inside  the  spore  coat  lays  the  outer  membrane.  Although  the  outer  membrane  is 
important  in  spore  formation,  in  the  dormant  spore  the  function  of  this  membrane  is  not 
clear  [33].  However,  like  the  spore  coat  the  outer  membrane  “has  no  notable  effect  of 
spore  resistance  to  heat”  [34]. 

liner  mi\  brane  Outer  msmbraie 


Figure  2.1:  Spore  Structure  [33] 

The  next  layer  in  the  spore  is  the  cortex  which  consists  of  peptidoglycan  (PG). 
The  cortex  reduces  the  core  water  content  which  is  essential  to  the  spore’s  resistance  to 
wet  heating.  The  germ  cell  wall  is  also  composed  of  PG.  During  germination,  the  germ 
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cell  wall  becomes  the  cell  wall  of  the  outgrowing  cell.  It  does  not  provide  any  heat 
resistance  for  the  spore  [12,  34], 

The  inner  membrane  separates  the  core  from  the  cortex  and  is  relatively 
impermeable  to  small  hydrophilic  and  hydrophobic  molecules  [33].  This  membrane 
provides  a  strong  permeability  barrier  to  protect  the  core  from  damage. 

In  the  center  of  the  spore  is  the  core.  The  core  consists  of  the  spore  DNA  and 
most  spore  enzymes.  Depending  on  the  species,  the  core’s  wet  weight  as  water  is  only 
25-55%.  Additionally,  the  amount  of  free  water  is  extremely  low  which  indicates  that 
molecular  mobility  within  the  core  is  almost  nonexistent.  The  core  water  content  plays  a 
key  role  in  spore  resistance  to  wet  heat  and  in  the  spore’s  enzymatic  dormancy.  The  spore 
is  more  resistance  to  wet  heat  with  low  core  water  content  [11].  The  core  water  content  is 
dependent  on  numerous  things  including  sporulation  condition,  temperature,  and  spore 
preparation  temperatures.  For  example,  the  higher  the  preparation  temperature,  the  lower 
the  core  water  content  [23].  In  high  humidity  conditions,  more  water  accumulates  in  the 
regions  external  to  the  core.  To  understand  the  chemical  reactions  taking  place  within  the 
spore,  reaction  kinetics  concepts  will  be  reviewed  in  the  next  section. 

Reaction  Kinetics 

Reaction  kinetics  will  be  used  throughout  to  describe  damage  reactions.  Reaction 
kinetics,  also  known  as  chemical  kinetics,  is  the  study  of  rates  of  chemical  processes. 
First-order  reactions  are  ones  that  are  directly  proportional  to  the  reactant  concentration. 
The  rate  of  the  reaction  is  directly  related  to  the  concentration  of  the  reactant  by  a 
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proportionality  constant,  k.  The  rate  constant  for  thermal  reactions  satisfies  an  Arrhenius 


equation,  that  is 


k{T)  =  Ae’^'^  (2.1) 

where  A  is  the  encounter  frequency,  Ea  is  the  activation  energy,  R  is  the  gas  constant,  and 
T  is  absolute  temperature.  Therefore,  the  rate  constants  increase  with  increasing 
temperature.  The  reaction  rates  are  influenced  by  the  physical  state  of  the  reactants,  the 
concentration  of  the  reactants,  the  temperature  at  which  the  reaction  occurs,  and  the 
presence  of  any  catalysts  in  the  reaction.  For  example,  later  in  Chapter  IV,  the  rate  of 
enzyme  damage  has  been  modeled  as 


dt 


-k,  (T)E  -  k^  T)E 


(2.2) 


where 


E  =  spore’s  enzyme  eoneentration 
[H2O]  =  average  water  concentration 

ki  =  rate  coefficient  associated  with  enzyme  breakdown  during  pyrolysis 
thermal  damage 

k2  =  rate  coefficient  associated  with  enzyme  breakdown  during  hydrolysis 
damage. 

The  units  of  enzyme  and  average  water  concentrations  are  in  molal  units.  Molal  is  the 
number  of  moles  of  solute  per  kilogram  of  solvent.  Activation  energies  for  hydrolysis 
reactions  are  generally  smaller  than  those  for  pyrolysis  reactions.  Therefore,  it  can  be 
seen  that,  since  water  plays  an  important  role  in  the  chemical  kinetics  of  thermal  damage, 
the  water  mobility  can  be  a  significant  influence  on  thermal  damage  within  the  spore. 
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Water  Mobility  Model 

Using  the  knowledge  of  the  role  of  each  layer  in  the  spore  with  respect  to  heat 
resistance.  Captain  Knight  developed  a  model  to  illustrate  the  mobility  of  water  within 
the  spore.  This  water  mobility  model  was  based  on  a  diffusion  equation  (Equation  2.3) 
where  only  the  spore’s  core  and  cortex  was  considered  in  a  two-layer  spore  model  [19], 
and  is  given  by 

—  =  /(h;)  +  V*Z)Vw  (2.3) 

dt 

where 

w  =  concentration  of  mobile  water  (molal) 
t  =  time  (sec) 

/(w)  =  production  of  water  (molal/sec) 

D  =  diffusion  coefficient  (cm^/sec). 

In  the  water  mobility  model,  since  the  core  initially  contains  an  extremely  low 
amount  of  absorbed  water,  it  was  assumed  that  the  concentration  of  absorbed  water  in  the 
core  is  zero.  Thus,  the  water  content  in  the  core  was  assumed  to  be  only  chemically 
bound  water.  Whereas,  the  cortex  contains  absorbed  and  chemically  bound  water 
molecules.  It  was  assumed  an  equal  amount  of  chemically  bound  water  is  available 
throughout  the  core  and  cortex. 

The  absorbed  and  chemically  bound  water  molecules  become  mobile  in  a  free 
molecular  state,  similar  to  a  gaseous  state,  when  the  spores  are  heated. 
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where 


[H2O]  I  =  initial  concentration  of  absorbed  water  molecules  (molal) 

[H2O]  b  =  initial  concentration  of  chemically  bound  water  molecules  (molal) 
[H2O]  g  =  initial  concentration  of  mobile  water  molecules  (molal). 

The  two  types  of  mobile  water  are  energized  at  different  temperatures.  These 
temperatures  will  be  estimated  in  the  Chapter  III  thermogravimetric  analysis. 

For  the  water  mobility  model,  the  initial  distribution  of  mobile  water  in  the  spore 
is  the  following, 

w(r,0)  =  [H,O\(r)  (2.4) 

where  r  is  the  radius. 

Since  water  is  consumed  through  reaction  with  hydrocarbon  compounds  within 
the  spore,  the  production  term  in  Equation  2.3  was  reduced  to  a  simple  first  order  reaction 
term: 

f(w)  =  -k,w  (2.5) 

where  kt  is  the  average  rate  for  all  possible  hydrolysis  reactions  in  the  spore.  The  rate 
coefficient  again  satisfies  the  Arrhenius  equation  given  in  Equation  2. 1 . 

In  the  water  mobility  model,  it  was  assumed  that  the  spore  is  spherically 
symmetric.  Therefore,  the  model  becomes  one  dimensional  with  respect  to  the  radius  of 
the  spore  and  time  (w(r,  t)).  The  diffusion  of  water  was  modeled  radially  from  the  spore’s 
coat  into  the  cortex  and  the  surrounding  environment.  At  the  boundary  where  the  radius 
is  equal  to  zero,  the  solution  must  be  finite,  that  is, 

I  w{0,t)  \<  00. 
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When  the  radius  is  equal  to  the  total  spore  radius  ( r  =  ),  the  flow  of  water  across  the 

outside  boundary  is  proportional  to  the  water  content  in  the  spore  at  the  outer  boundary 
and  ambient  humidity.  This  boundary  condition  was  modeled  by, 

-h  ■  DVw(ro  ,t)  =  h{  w{r^ ,  t)  -  [H^O\ }  (2.6) 

where 

h  =  outward  normal  vector 
[H2O]  a  =  water  concentration  outside  spore  (molal) 
h  =  surface  flow  velocity. 

For  spherical  geometry,  h  became  the  outward  radial  direction.  Therefore,  Equation  2.4 
became 


-D^{r,,t)  =  h{w{r„t)-{H,0\}. 
dr 


(2.7) 


Since  the  diffusivity,  D,  was  assumed  to  remain  constant  across  the  radius  of  the 
spore,  the  diffusion  term  in  Equation  2.3  became 


V  DS/w  =  D 


d^w  2  dw 
dr^  r  dr 


(2.8) 


In  order  to  find  the  analytical  solution  to  the  water  mobility  model,  dimensionless 
variables  were  chosen.  Refer  to  Ref.  19  for  additional  details  pertaining  to  these 
variables.  Introducing  these  variables  gives  the  following  expression, 

—  =  -^v  +  ^  — (x"— )  0<x<l,r>0  (2.9) 

dr  dx  dx 


with  boundary  conditions 


—  (0,r)  =  0  r>0  (2.10) 

dx 
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r>0  (2.11) 


-^(l,r)  =  5(v(l,r)-7) 
ox 


The  initial  condition 


v(x,0)  =  Qrj(x)/j(r*)  +  Qr2W/2(7’*)  0<x<l  (2.12) 

The  dimensionless  constants  are 


R 


P  = 


D 


and 


iH,Ol 


with  functions 


a^(x)  = 


[H,Ol(r,x) 


and 


ajix)  = 


[H,Ol(r,x) 


Using  separation  of  variables,  the  homogenous  problem  given  in  Equation  (2.7) 
and  the  homogenous  boundary  conditions  became 


/^'(T)  +  (/?  +  //)/^(^)  =  0 

(2.13) 

[x^g'(x)]’-h//x^g(x)  =  0 

(2.14) 

g'(0)  =  0 

(2.15) 

-g'(l)  =  5g(l) 

(2.16) 

where  ju  is  a  separation  parameter.  The  eigenvalues  and  eigenfunctions  were  found  by 
solving  Equations  2.14,  2.15,  and  2.16  with  respect  to  the  homogenous  boundary 
conditions.  Further  details  are  outlined  in  Ref.  19  page  3-14  and  3-15. 

The  solution  to  the  homogenous  system  is  the  following. 


00  00  2 

v(x.T)  =  j;,Fjr)e-''-’g,(x)  +  ry{l,g.)^li-e-'-']g„U)-  (2.17) 

«=1  n=l 

The  first  series  represents  the  influence  of  the  initial  water  distribution  while  the  second 
series  represents  the  influence  of  the  wet  or  dry  heating  environment.  Since  the  DNA  is 
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only  present  in  the  spore  core,  all  hydrolysis  reactions  with  DNA  will  occur  in  the  core. 
The  average  total  mobile  water  concentration  in  the  core  was  given  by 


{H^Ol{t)  =  -^lw{r,tydr.  (2.18) 

Substituting  Equation  2.17  and  dimensionless  variables  into  Equation  2.18  yields 

O  f  CO  CO  1 


n-l 


Physically,  Equation  2.19  states  that  for  short  exposure  times,  the  second  term 
approaches  zero.  This  implies  that  the  first  term  dictates  the  release,  diffusion,  and 
reaction  of  water  in  the  core.  During  long  exposure  times,  the  effect  of  the  water 
concentration  outside  the  spore  becomes  relevant  [19].  Eor  additional  details  about  the 
water  mobility  model  see  Ref.  19  Appendix  B. 

Captain  Knight’s  water  mobility  model  will  be  used  to  characterize  the  spore’s 
chances  of  survival  during  the  inactivation  process.  The  spore’s  survival  is  dependent  on 
hydrolysis  and  pyrolysis  reactions  occurring  within  the  spore  during  the  thermal  heating. 
Chapter  III  Thermogravimetric  Analysis  will  focus  on  the  concentration  of  water  within 
the  spore  and  the  temperatures  at  which  this  water  will  be  available  for  hydrolysis 
reactions. 
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III.  Thermogravimetric  Analysis 


Introduction 

Two  types  of  water  are  contained  within  the  spore,  absorbed  and  chemically 
bound  water.  Absorbed  water,  weakly  associated  with  cell  materials,  is  located  in  the 
spore’s  cortex  while  chemically  bound  water  is  located  in  the  spore’s  core  and  cortex. 
During  the  inactivation  of  the  spore,  the  absorbed  and  chemically  bound  water  molecules 
become  mobile.  The  weakly  bound  absorbed  water  becomes  energized  and  mobile  at 
lower  temperature  than  strongly  held  chemically  bound  water.  Captain  Knight  assumed 
these  temperatures  were  around  100  °C  for  absorbed  water  and  300  °C  for  chemically 
bound  water  [19].  In  order  to  confirm  these  assumptions  and  approximate  the  percentage 
of  water  in  the  Bacillus  spore,  thermogravimetric  analysis  studies  were  conducted  by  Dr. 
Felker  [10].  Dr.  Felker  used  Bacillus  thuringiensis  when  conducting  this  analysis.  The 
percentage  of  water  in  the  spores  is  important  during  the  inactivation  process;  the 
diffusion  of  water  allows  hydrolysis  reactions  to  occur  within  the  spore.  These  reactions 
are  important  in  the  inactivation  of  the  spore  because  they  lead  to  deploymerization 
reactions  of  the  spore’s  DNA  and  enzymes.  Damage  to  the  DNA  and  to  enzymes 
influence  the  DNA  repair  process  and  the  spore’s  ability  to  produce  outgrowth  as  well  as 
the  repair  time  to  produce  outgrowth. 

General  Discussion 

Thermogravimetric  analysis  (TGA)  is  a  thermal  analysis  technique  that  measures 
the  weight  change  in  a  material  as  a  function  of  temperature  and  time  under  a  controlled 
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atmosphere  [38].  Thus,  TGA  can  directly  measure  the  weight  loss  of  a  material  with 
temperature  or  time,  as  such  loss  is  due  to  dehydration  or  decomposition  as  the 
temperatures  ranges  from  ambient  to  1000°C.  TGA  can  be  used  in  multiple  applications, 
such  as,  thermal  stability  or  degradation  investigation  of  organic/inorganic  materials  and 
inert/oxidative  atmospheres,  determining  of  organic/inorganic  content  of  mixtures,  phase 
transition  measurement,  reaction  kinetics  with  reactive  gases,  and  pyrolysis  kinetics  [38]. 

The  TGA  of  solids  relies  on  a  high  degree  of  precision  in  three  measurements: 
weight,  temperature,  and  temperature  change.  The  analysis  is  carried  out  by  raising  the 
temperature  gradually  and  plotting  percentage  of  weight  change  against  temperature  [28, 
38].  The  thermogravimetric  analyzer  consists  of  a  high  precision  balance  with  a  platinum 
pan  loaded  with  a  sample.  The  platinum  pan  is  rigidly  attached  to  the  end  of  the  balance 
arm  in  the  horizontal  configuration,  and  there  is  no  added  mass  due  to  a  suspension 
harness.  According  to  Menczel  et  al,  “suspension  below  the  balance  may  be  less  rigid  and 
does  not  recognize  the  compensating  weight  and  thus  detracts  less  from  the  available 
range”  [24].  The  sample  is  placed  in  a  small  electrically  heated  oven  with  a  thermocouple 
to  accurately  measure  the  temperature.  A  computer  is  used  to  control  the  instrument. 

TGA  can  provide  a  quantitative  measurement  of  the  weight  change  but  cannot 
indicate  the  nature  of  the  material  lost.  Therefore,  evolved  gas  analysis  (EGA) 
measurements  can  be  done.  EGA  gives  the  nature  of  the  changes  in  the  gas  phase.  It  can 
be  qualitative  or  quantitative  measurements.  Mass  spectroscopy  (MS)  or  Eourier 
transform  infrared  spectroscopy  (ETIR)  can  be  applied  to  the  EGA  [24].  Eor  our  analysis 
Dr.  Eelker  used  ETIR.  The  ETIR  system  allowed  Dr.  Eelker  to  confirm  that  the  major 
volatile  associated  with  the  mass  change  was  indeed  water. 
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Thermogravimetric  Analysis  Data 

Raw  TGA  Data. 

There  are  two  types  of  TGA,  differential  and  normal.  The  normal  TGA  method 
was  used  for  our  study.  The  sample  can  be  in  various  atmospheres  including  oxidizing 
and  non-oxidizing.  An  oxidizing  atmosphere  is  one  in  which  oxygen  is  present  whereas 
in  a  non-oxidizing  atmosphere  oxygen  is  not  present.  For  this  study,  Bacillus 
thuringiensis  (B.t)  was  focused  on.  The  B.t  sample  was  in  a  quasi-wet  state  and  slowly 
dried.  Therefore,  the  study  has  a  steady  state  measurement  of  all  excess  water.  Although, 
there  is  water  located  on  the  outside  of  the  sample  and  within  the  spores,  this  study  is 
only  concerned  with  the  water  within  the  spores.  As  the  sample  is  heated,  several  things 
can  happen.  The  sample  can  stay  the  same,  lose  mass,  or  gain  mass.  Also,  as  the  sample 
is  heated,  more  than  one  process  is  occurring.  There  is  an  early  release  of  water  and  at 
later  temperatures  combustion  is  released  from  the  sample. 

The  TGA  raw  data  obtained  from  using  B.t  can  be  seen  in  Figure  3.1.  The  blue 
curve  represents  the  weight  change  in  the  sample  as  a  function  of  time  (seconds).  The 
green  curve  shows  the  increase  in  temperature  with  time.  The  temperature  increases 
linearly  at  a  rate  of  0.2  °C  per  second.  Also,  there  are  distinct  slope  variations  in  the  mass 
as  time  increases. 
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For  use  in  this  study,  the  TGA  data  was  represented  by  mass  change  as  a  function 
of  temperature.  This  is  the  blue  curve  depicted  in  Figure  3.2.  This  curve  depicts  steep 
slopes  between  55°C  and  125°C  and  between  200°C  and  350°C.  These  steep  slopes 
indicate  a  rapid  loss  of  mass  in  these  temperature  ranges.  This  mass  loss  is  related  to  the 
release  of  water  from  the  sample. 
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Modeling  Data. 

In  order  to  estimate  the  amount  of  absorbed  and  chemically  bound  water  released 
from  spores  at  selected  temperatures,  a  model  was  built  to  provide  the  best  fit  to  the  TGA 
experimental  data  (blue  curve  in  Figure  3.2).  Gaussian  density  functions  were  chosen  to 
represent  the  distributions  of  absorbed  water,  chemically  bound  water,  and  other 
chemical  processes.  Each  Gaussian  density  function  is  centered  about  their  respective 
mean  release  temperature,  Ti,  T2,  and  T3  with  associated  standard  deviations  cti,  CT2,  and 

^3. 
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.  Ailzk 

f,{T)  = - 


/  =  1,2,3. 


The  total  mass,  G(T  ),  released  at  any  given  temperature,  T  ,  is  determined  by 
integrating  the  weighted  sum  of  the  Gaussian  density  functions,  that  is. 


G(T*)=  f  Wj/j(r)  +  W2/2(r)  +  W3/3(r)  dr 

J-oo 


=  w.F.iT*)  +  w.F.iT*)  +  w.F.iT*) 


(3.1) 


where 


F,(r)= 


1 

2 


l  +  erf 


^  cr,V2  J_’ 


/  =  1,2,3. 


(3.2) 


Here  wi  represents  the  total  amount  of  absorbed  water  available  for  release,  W2  represents 
the  total  amount  of  chemically  bound  water  available  for  release,  W3  represents  the  total 
amount  of  other  chemical  processes  available  for  release,  and  erf  is  the  error  function. 

The  TGA  model  is  then  expressed  by  subtracting  the  total  mass  released  at  any 
temperature,  G(T*),  from  the  total  mass  of  the  sample,  TM,  i.e. 

TGAf  (T*  )  =  TM-  G{T* )  (3.3) 


where 

TM  =  Wj  +W2  +W3. 

Furthermore,  how  much  water  is  released  at  any  temperature  from  the  spore, 
[H20]s(T*),  can  be  determined  by  subtracting  the  total  available  water,  [H20]a,  released 
at  any  temperature  from  the  total  mass  of  the  sample. 
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{H^Ol{r)  =  TM  -{H^OUT*) 


(3.4) 


where 


{H,0\=w,F,{r)  +  w,F,{r). 


Results. 

Using  the  TGA  model,  the  appropriate  values  that  give  the  best  fit  approximation 
to  the  TGA  experimental  data  curve.  These  values  are  outlined  below  in  Table  3.1. 


Table  3.1:  Parameters  for  TGA  Curve  Fit 


Weiahted 

Mean 

Standard 

Coefficient 

Temneratnre 

Deviation 

(Wi) 

(Ti) 

Absorbed  Water 

5 

105 

60 

Chemically  Bound  Water 

38 

290 

70 

Other  Chemical  Processes 

50 

815 

300 

The  three  weighted  Gaussian  density  functions  are  illustrated  in  the  Figure  3.3 
using  the  values  from  Table  3.1.  The  majority  of  absorbed  water  is  mostly  likely  released 
in  a  neighborhood  of  105  “C  until  all  the  absorbed  water  is  active.  Then  the  chemically 
bound  water  most  likely  begins  its  release  to  a  peak  of  290  °C.  The  other  chemical 
processes  which  reduce  the  mass  are  mainly  released  to  a  peak  of  815  °C. 
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Temperature  (Degrees  C) 

Figure  3.3:  Weighted  Gaussian  Density  Functions 

Figure  3.4  shows  the  results  of  our  best  fit  approximation  to  the  TGA  data  using 
the  values  in  Table  3.1.  This  figure  depicts  the  release  of  mass  from  the  sample  as 
temperature  increases.  The  blue  curve  represents  the  true  TGA  experimental  data,  the 
solid  red  line  represents  our  best  fit  model  to  the  TGA  data,  and  the  dashed  red  line 
represents  the  amount  of  mass  loss  due  to  water  only. 

From  probability  theory,  the  empirical  rule  states  that  for  a  Gaussian  distribution 

[12], 

68%  of  the  measurements  lie  within  fx±a 
95%  of  the  measurements  lie  within  jj.  ±  2cr 
97.7%  of  the  measurements  lie  within  //±3f7 
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where  )a,  represents  the  mean  and  a  represents  the  standard  deviation.  Therefore,  since  our 
best  fit  model  gives  a  close  approximation  to  the  TGA  experimental  data,  68%  of 
absorbed  water  is  released  approximately  between  45°C  and  165°C  while  68%  of  the 
chemically  bound  water  is  released  approximately  between  220°C  and  ShO^C.  Also,  our 
results  support  Captain  Knight,  assumptions  that  the  absorbed  water  becomes  available 
around  100°C  and  chemically  bound  water  becomes  available  around  300°C. 


Figure  3.4:  TGA  Model  Fit 

Figure  3.5  shows  the  percent  mass  loss  as  temperature  increases.  Approximately 
2%  of  the  sample’s  weight  is  loss  around  105°C  and  25%  is  loss  around  290°C. 
Furthermore,  by  430°C  all  of  the  sample’s  water  is  released  and  that  approximately  40- 
45%  of  the  sample’s  weight  is  water.  The  percentage  of  water  within  the  spore  is 
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important  for  hydrolysis  reactions  to  occur.  The  higher  percentage  of  water  implies 
increase  hydrolysis  reactions  can  occur.  During  the  heating  process  of  the  Bacillus 
spores,  hydrolysis  reactions  are  an  important  aspect  to  inactivating  the  spores. 


Figure  3.5:  TGA  Model  Fit  (Percent  Mass  Loss) 
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IV.  Spore  Damage  and  Repair  Mechanisms 


Introduction 

In  Chapter  III,  the  study  observed  that  approximately  40-45%  of  the  spore’s 
weight  is  water.  Water  plays  an  important  role  in  the  inactivation  of  the  spores.  The 
transport  of  water  throughout  the  spore  allows  hydrolysis  reactions  to  occur  in  the  core. 
The  hydrolysis  reactions  can  cause  damage  to  the  spore  DNA  and  enzymes  which  leads 
to  the  spore  death  if  the  DNA  is  not  repaired.  In  the  spore  population,  it  is  assume  that 
each  spore  contains  a  full  complement  of  DNA  and  a  concentration  of  enzymes  adequate 
for  germination.  It  has  been  assumed  that  the  damage  repair  capability  for  a  spore  in  the 
spore  population  is  a  sample  from  a  certain  initial  enzyme  fitness  distribution.  When  the 
enzyme  damage  accumulates,  this  fitness  distribution  is  degraded.  Once  degradation  is 
beyond  some  critical  level,  the  spore  cannot  repair  the  damaged  DNA  and  therefore, 
cannot  germinate  or  produce  outgrowth.  Probability  of  kill  models  based  on  water 
mobility,  hydrolysis,  DNA  information  content,  and  the  spore’s  enzyme  fitness  levels 
have  been  created  to  further  explore  the  spore’s  chance  of  survival. 

Spore  Damage  Model 

Hydrolysis  &  Pyrolysis. 

Hydrolysis  is  the  reaction  of  a  molecule  with  water.  At  low  temperatures, 
hydrolysis  reactions  can  occur  because  of  their  low  activation  energy.  Hydrolysis  can  be 
catalyzed  by  hydrogen  ions  in  low  pH  or  by  hydroxide  ions  in  high  pH,  where  the  pH  of 
an  aqueous  solution  is  the  logarithm  of  the  reciprocal  of  the  concentration  of  hydrogen 
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ions  in  the  solution.  A  solution  with  a  pH  below  seven  is  acidic  and  a  pH  above  seven  is 
basic;  neutral  solution  has  pH  equal  to  seven.  So  hydrolysis  rates  increase  in  strongly 
basic  or  strongly  acidic  solution. 

(Note  that:  [//  ]  +  k(jjg^y.gjj[H20]  +  ])[32] 

Bacteria  can  survive  in  an  environment  with  pH  measurements  between  five  and 
eight,  so  acid  hydrolysis  is  somewhat  less  damaging  than  basic  hydrolysis  paralleling  the 
pH  dependence  of  hydrolysis  rates  for  most  organic  molecules  that  are  susceptible  to 

hydrolysis  for  which  k(acid  catalysis)  <  k(base  catalysis). 

Autoionization  is  a  reaction  in  which  water  molecules  dissociate  into  hydrogen 
cations  (H"^)  and  hydroxide  anions  (OH  )  in  an  equilibrium  process.  That  is, 

+OH-. 

The  equilibrium  constant  (K)  for  this  dissociation  of  water  is  defined  as 

The  pH  in  a  solution  is  related  to  its  acidity  or  basicity.  pH  is  defined  as 

pH  =  -log[/f"] 

where  [H"^]  is  in  molarity  (M).  The  concentration  of  hydroxide  ions  in  the  solution,  given 
by 

pOH=-log[6>H-] 

can  be  derived  from  pH  using  the  equilibrium  constant  (K).  Define  a  new  constant,  Kw, 
such  that 
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Then 


Therefore, 


^  [0//-]= 


[//"] 


-log[0//-]  =  -log-^, 

=  -\ogK^+log[H^l 


^  pOH  =  -  log  -  pH. 

The  change  in  pH  within  the  spore  can  cause  damage  to  the  spore  DNA  by 
hydrolysis.  For  example,  as  pH  approaches  9,  the  DNA  double  helix  begins  to  dissociate 
into  its  component  single  strands.  This  is  because  the  hydroxide  ions  can  react  with  bases 
in  DNA  base  pairs  to  remove  certain  protons  enhancing  the  rate  of  reaction  of  the  bases 
with  water.  Acids  can  accelerate  the  rate  of  reaction  of  DNA  bases  with  water  by 
donating  a  proton  to  the  bases.  Therefore,  certain  chemicals  that  alter  spore  pH  may  be 
effective  to  cause  DNA  damage  which  can  lead  to  spore  death. 

Pyrolysis  is  the  thermal  decomposition  of  organic  materials  within  the  spores. 
These  damage  reactions  are  mediated  by  radical  formation,  abstraction  and 
recombination  reactions.  Pyrolysis  usually  occurs  at  much  higher  temperatures  than  do 
hydrolysis  reactions  except  in  the  extreme  absence  of  water. 

Types  of  Damages  to  DNA  and  enzymes. 

DNA  and  enzymes  are  biopolymers  having  very  different  structures  and  very 
different  functions.  DNA  is  the  information  storage  medium  that  encodes  the  recipe  for 
making  all  cellular  polymers,  including  the  enzymes.  Enzymes  are  active  proteins  which 
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catalyze  specific  chemical  reactions  necessary  to  build  biopolymers  and  disassemble 
them  to  generate  energy. 

DNA  is  a  linear  polymer  made  up  of  four  different  types  of  monomers.  It  has  a 
fixed  backbone  which  is  built  of  repeating  sugar-phosphate  units.  Each  sugar  is 
connected  to  one  of  four  possible  bases:  adenine  (A),  cytosine  (C),  guanine  (G),  and 
thymine  (T).  In  1953,  James  Watson  and  Francis  Crick  correctly  proposed  that  the  DNA 
molecule  consist  of  two  stands  bound  together  in  a  double  helix  through  specific 
hydrogen  bonding  affinities  of  complementary  bases  for  each  other.  The  double  helix 
structure  is  arranged  such  that  the  sugar-phosphate  backbone  lies  on  the  outside  and  the 
bases  on  the  inside.  The  bases  form  specific  base  pairing  held  together  by  hydrogen 
bonds:  adenine  (A)  pairs  with  thymine  (T)  and  guanine  (G)  pairs  with  cytosine  (C).  The 
sequence  of  bases  along  the  DNA  strands  is  the  manner  that  genetic  information  is  stored. 
The  DNA  sequence  also  determines  the  sequences  of  the  ribonucleic  acid  (RNA)  and 
protein  molecules.  If  the  DNA  double  helix  is  separated  into  two  single  strands,  each 
strand  can  act  as  a  template  for  replication.  For  perfect  encoding  of  information,  the  loss 
of  a  single  base  in  one  DNA  strand  can  be  recognized  by  reading  the  complementary  base 
and  so  the  total  information  in  the  double  helix  is  not  degraded  if  the  loss  is  recognized 
and  repaired. 

During  the  thermal  inactivation  of  the  spores,  there  are  numerous  types  of  reactive 
processes  that  can  lead  to  DNA  damage.  Some  include  chemical  modification  such  as 
alkylation  or  oxidation  and  base  removal  [35].  However,  not  all  regions  of  DNA  within 
the  spore  are  important  for  spore  survival  and  replication.  Therefore,  only  damage  to  the 
critical  or  essential  DNA  necessary  for  germination  and  outgrowth  is  considered. 
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Proteins  are  linear  polymers.  Proteins  are  built  from  20  amino  acids  and  fold  up 
into  elaborate  three-dimensional  structures  determined  by  their  amino  acid  sequences. 
Enzymes  are  proteins  folded  into  three-dimensional  structures  that  can  catalyze  specific 
chemical  reactions.  DNA  genes  encode  the  sequences  of  proteins.  The  gene  is  the 
fundamental  unit  of  hereditary  information  encoding  a  series  of  3 -bit  letters  or  codons.  A 
codon,  a  set  of  three  bases  along  the  DNA  strand,  determines  the  identity  of  one  amino 
acid  within  the  protein  sequence.  One  incorrect  amino  acid  can  dramatically  alter  the 
function  of  the  protein.  The  three  bits  of  ordered  information  allow  64  possible  encoding 
combinations.  The  capability  of  the  three-bit  letter  is  sufficient  to  encode  for  the  20 
amino  acids.  Therefore,  there  is  redundancy  in  the  code.  As  a  result,  removing  or  altering 
the  third  base  pair  in  the  sequence  in  a  codon  may  or  may  not  destroy  the  information 
content.  Whereas,  removing  the  first  base  in  a  codon  most  likely  will  cause  problems 
[19].  Some  codon  letters  are  inactive,  so  destroying  that  information  will  cause  no 
degradation  in  the  useful  DNA  gene  information  or  proteins  function  encoded  by  that 
gene. 

As  previously  stated,  spore  DNA  and  enzyme  damage  can  occur  in  many  ways. 
According  to  Setlow  “chemicals  such  as  nitrous  acid  or  formaldehyde  can  kill  the  spores 
by  DNA  damage,  while  others,  such  as  oxidizing  agents,  appear  to  damage  the  spore’s 
inner  membrane  so  that  this  membrane  ruptures  upon  spore  germination  and  outgrowth” 
[34].  Water  can  react  with  DNA  and  proteins  resulting  in  deploymerization. 
Deploymerization  is  breaking  a  large  strand  of  polymers  into  two  smaller  ones.  Water 
more  frequently  reacts  with  DNA  causing  deamination  and  depurination  damage. 
Deamination  is  the  removal  of  amino  group  from  the  DNA’s  base  pair.  This  type  of 
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damage  can  prevent  the  DNA  from  replicating  properly.  The  major  type  of  deamination 
reaction  converts  cytosine  (C)  to  an  altered  DNA  base,  uracil  (U),  but  deamination  occurs 
on  other  bases  as  well. 

Denaturation  is  the  reaction  in  which  the  proteins  lose  their  structure.  The  proteins 
lose  their  structure  by  application  of  some  external  stress,  such  as  heat  or  a  strong  acid  or 
base.  Following  denaturation,  however,  some  proteins  will  return  to  their  native  structure 
under  proper  conditions.  Although,  under  extreme  conditions  usually  cause  irreversible 
change  and  the  protein  can  no  longer  perform  its  function  once  it  has  been  denatured. 

Depurination  is  an  alteration  of  DNA  in  which  the  purine  base  (adenine  or 
guanine)  is  removed  from  the  DNA  chain.  Depurination  damage  is  most  common  since 
the  purine  bond  is  especially  susceptible  to  hydrolysis. 

When  DNA  is  missing  one  of  its  base  pairs,  the  DNA  strand  cannot  match  up  with 
its  other  half  of  the  helix  chain.  Removal  of  a  DNA  base  can  have  a  range  of  effects 
depending  on  the  information  content  of  the  codon.  For  this  research,  only  damage  which 
reduces  the  information  content  in  critical  or  essential  codon,  which  contain  information 
necessary  to  replicate  the  spore  during  germination.  Also,  it  has  been  assumed  that 
hydrolysis  reactions  are  equally  likely  to  occur  at  all  protein  sites  and  that  the  critical 
threshold  damage  to  DNA  caused  by  depurination  is  equally  likely  for  all  base  pair  sites 
[19]. 

DNA  Repair  Mechanisms. 

DNA  repair  is  a  complex  system  of  interacting  biomolecular  processes.  The  DNA 
repair  pathways  have  evolved  but  are  partially  conserved  in  all  life,  from  bacteria  to 
mammals  [36].  Each  repair  pathway  is  directed  to  a  specific  type  of  damage,  and  several 
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pathways  can  target  a  given  type  of  damage.  Major  DNA  repair  pathways  include  base 
excision  repair  (BER),  nucleotide  excision  repair  (NER),  mismatch  repair  (MMR), 
homologous  recombinational  repair  (HR),  and  non-homologous  end  joining  (NHEJ). 
Each  repair  pathway  requires  a  number  of  repair  enzymes  for  the  process. 

BER  is  the  main  mechanism  for  removing  spontaneous  DNA  lesions  that  are 
caused  by  hydrolysis  or  oxidation.  This  repair  process  removes  damaged  bases  from 
DNA  through  a  multistep  process  in  which  the  repair  enzymes  work  independent  of  each 
other  [37].  An  example  of  BER  is  the  repair  of  uracil-containing  DNA.  Recall  that  uracil 
can  be  formed  in  DNA  by  the  deamination  of  cytosine. 

NER  removes  bulky  adducts  in  DNA  or  chemical  bonds  between  bases  of  DNA. 
The  mechanisms  of  NER  are  more  complex  than  that  of  BER,  but  the  basic  principles  are 
similar  [37].  MMR  provides  protection  against  occasional  rare  errors  that  occur  during 
DNA  replication  [8].  HR  and  NHEJ  double-stand  breaks.  HR  uses  a  homologous  DNA 
template,  is  highly  accurate,  whereas  NHEJ  rejoins  the  broken  ends  without  a  template 
and  is  often  accompanied  by  loss  of  some  nucleotides  [8]. 

As  mentioned  above,  deamination  and  depurination  is  the  reaction  of  water  with 
the  base  pairs  of  DNA.  Therefore,  most  likely  during  the  thermal  inactivation  of  the 
spores,  the  base  pairs  are  damaged.  Since  BER  is  the  main  repair  system  for  this  type  of 
damage,  it  will  be  assumed  during  germination  that  BER  will  be  used  to  repair  most  of 
the  damaged  DNA,  dominating  the  overall  rate  of  repair.  Eurther  discussion  will  be  given 
in  Chapter  V. 
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Damage  Model. 

This  research  assumes  that  the  spore  starts  out  with  100%  healthy  DNA  and  this  is 
degraded  by  exposure  to  heat  due  to  hydrolysis  and  pyrolysis  reactions.  The  rate  of 
critical  DNA  damage  can  be  modeled  by  the  following  rate  equation, 

=  k,{T)I^  (4.1) 

dt, 

where represents  the  critical  DNA,  ^jis  the  rate  coefficient  associated  with  critical 

DNA  breakdown  during  pyrolysis,  ^2  is  the  rate  coefficient  associated  with  critical  DNA 

breakdown  during  hydrolysis,  te  is  the  exposure  time,  T  is  temperature,  and  [H20]c  is  the 
average  core  water  concentration  given  in  Equation  2.19.  The  rate  coefficients  again 
satisfy  the  Arrhenius  equation  in  Equation  2.1. 

A  simple  solution  to  the  rate  equation  is  the  following  where  the  critical  DNA 
diminishes  exponentially, 

IoiO  =  IMT,0  (4.2) 

where  =  the  fraction  of  DNA  remaining  after  some  exposure 

time  and  temperature,  H{t^)=  ,  and  Ic  represents  the  initial  state  of  critical 

DNA. 

As  the  spore  is  heated,  there  is  also  damage  occurring  to  the  spore’s  enzymes. 
Since  the  enzymes  are  crucial  for  DNA  repair,  a  damaged  spore  will  only  be  able  to 
germinate  using  the  enzymes  that  survive  the  heating  process.  The  damage  to  the  DNA 
must  be  repaired  before  a  viable  cell  is  produced  during  outgrowth  and  there  has  to  be  a 
suitable  amount  of  repair  enzymes  viable  to  accomplish  this. 
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Let£'^(?^)  be  the  enzyme  concentration  for  the  spore  population  at  exposure, 


The  rate  of  enzyme  damage  can  be  modeled  by 

dF 

(4.3) 

dt^ 

where  ki  is  the  rate  coefficient  associated  with  enzyme  breakdown  during  pyrolysis  and 
k2  is  the  rate  coefficient  associated  with  enzyme  breakdown  during  hydrolysis.  These  rate 
coefficients  may  not  be  the  same  rate  coefficients  associated  with  DNA  damage.  After 
integration.  Equation  4.3  becomes 

=  E,f{T,Q  (4.4) 


where,  for  simplicity, 

=  (4.5) 

Equation  4.5  represents  the  fraction  of  remaining  enzymes  after  a  given  exposure  time 
and  temperature. 


Fitness  Distribution  Models 

The  spores’  fitness  distribution  has  been  modeled  with  a  Gaussian  distribution 
representation  [19].  In  this  section,  this  approach  will  be  reviewed  and  extended  the 
fitness  model  to  a  beta  distribution  representation. 

Gaussian  Model. 

The  Gaussian  (Normal)  distribution  is  the  most  widely  used  continuous 
probability  distribution  [12].  The  Gaussian  probability  density  function  has  the  following 
representation  with  mean,  ju  and  variance,  , 
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f(x)  =  — 7^^  5  -oo<x<oo,  -co<ju<co,  (7>0.  (4.6) 

errin' 

The  Gaussian  cumulative  density  function  can  be  found  by  integrating  equation  4.5, 

{‘’e  dx.  (4.7) 

J-oo 

A  very  important  contributor  to  spore  fitness  is  the  store  of  enzymes  in  the  spore 
that  can  serve  to  reconstitute  an  active  cell.  Prior  to  any  heat  treatment,  the  initial  enzyme 
content.  Eg,  will  be  dependent  on  the  initial  state  of  the  spore  population.  This  initial  state 
of  enzymes  will  depend  on  numerous  things.  For  example,  spore  preparation  or 
sporulation  will  have  an  important  role  on  the  initial  state.  According  to  Melly  et  al,  spore 
preparation  temperature  affects  the  core  water  content  which  in  turn  influences  the 
spore’s  heat  resistance  [23].  Therefore,  was  drawn  from  a  Gaussian  fitness  distribution 

with  mean,  and  variance,  cjj,^ , 

1  2cr  ^ 

f(EQ)  = - 7^^  ^  ’  -00<£'q<00,  -co</j.^<co,  >0. 

CT 77  ^ '^TT 

The  illustration  of  a  notional  enzyme  fitness  density  function  prior  to  heat  is  shown  in 
Figure  4.1. 

In  order  to  repair  damaged  DNA,  a  certain  amount  of  viable  enzymes  are  needed. 
However,  during  the  inactivation  of  the  spore,  damage  is  also  done  to  the  enzymes. 
Therefore,  there  must  be  a  point  at  which  there  will  be  insufficient  concentration  of 
enzymes  to  repair  the  damaged  DNA.  This  location  is  the  critical  kill  threshold,  Ek. 
Above  this  threshold,  the  spore  contains  the  necessary  enzymes  needed  for  repair. 


P(x<Xq)  =  <l)(x)  =- 
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Figure  4.1:  Enzyme  Fitness  Prior  to  Any  Heating 

As  thermal  exposure  time  progresses,  the  evolution  of  the  distribution  of 
£'^(t^)will  approaeh  the  kill  threshold.  Since  Fq  is  a  random  variable,  )  is  also  a 

random  variable  drawn  from  a  distribution  with  mean,  and  variance,  cr^  ^ . 

From  probability  theory,  the  expectation  of  a  continuous  variable  X  is 

E(X)=  r  xf(x)dx,  (4.8) 

J-oo 

where  f(x)  is  the  probability  density  function  for  X.  The  mean  and  variance  of  a 

continuous  distribution  is  defined  in  terms  of  expectations,  that  is 

M  =  E[X] 
a^=E[(X-^f] 
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Expectation  is  a  linear  operator,  that  is 

E{cX]  =  cEiX] 

+ ...  +  X  J  =  £[XJ  +  £[X  J  + ...  +  E[X  J 

where  c  is  a  constant  and  Xi,  X2,. . .,  Xn  are  continuous  random  variables  [12]. 
Using  these  expectation  properties,  the  mean  of  E^{t^)  is  the  following 


=E{EM^)} 

=  EiEJ,{T,Q-\ 

and  variance  of  £'^(t^)is  the  following 


Using  Equation  4.8, 


=D^« 

Eq 


dEn 


dEr, 


Using  integration  by  parts  method,  ^udv  =  u-\- ^vdu  ,  let 


u  =  ^  °  dv  =  (E, )e  dE„ 

(7p  yl'2.7r 


then 
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du  = 


1 


dEr, 


<7p  ^27r 


v  =  -a„  e 


2„ 


Thus, 


Op  ^/^ 


2~  2(a^n)' 


(Jr  e 


+  - 


(Jr  e 


dEr, 


=  0  +  (o-£  f  r - 


(Eo-MEpf 
'  ^(^Epf 


dE„ 


=  (o-E.f- 


Therefore, 


<yl=ifEiT,K)f  E[{E,-pip^f] 


Figure  4.2  depicts  the  evolution  of  the  distribution  of  Er(te)  as  exposure  time 
increases.  The  Gaussian  function  outlined  in  black  represents  the  initial  distribution  of  the 
enzyme  content  with  mean,  and  variance,  c^Ep  ■  The  Gaussian  function  outlined  in 
blue  represents  the  distribution  of  the  enzyme  content  after  some  exposure  time,  t  =  z^ 
with  mean, (Tj)  and  variance,  ^(Tj)  .  The  Gaussian  function  outlined  in  red 
represents  the  distribution  of  the  enzyme  content  after  additional  heat 
exposure,  t  =  Tj  >  Tj  with  mean,//^  (tj)  and  variance,  cr^  ^(72) .  As  the  exposure  time 
increases,  the  mean  of  the  distribution  decreases  as  does  the  variance. 
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Figure  4.2:  Enzyme  Fitness  Evolution 

The  probability  of  kill,  Pk,  for  the  spore  population  is  the  cumulative  probability 
up  to  the  kill  threshold,  ^  ) ,  that  is 

=  P,{EJ,{T,0<E,) 

^  E  ^ 

Eo<  ' 


V 


fEiT,0 


J 


Using  equation  4.7, 


Eh 

ffE(T,te 

J-oo 


(Eo-MEpf 


dE^ 


=  \\l+erf{Q,yi 
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where 


Beta  Model. 

The  Gaussian  distribution  is  often  used  because  the  random  variables  can  be 
standardized,  but  this  assumes  the  random  variables  are  between  negative  infinity  and 
positive  infinity.  However,  since  the  maximum  amount  of  enzymes,  Emax,  within  any 
spore  cannot  be  infinite  nor  can  the  minimum  amount  be  negative,  the  Gaussian 
representation  may  not  give  an  adequate  physical  representation.  A  probability 
distribution  bounded  by  [0,  Emax]  would  provide  a  better  estimation  of  the  fitness 
probability.  Therefore,  the  beta  distribution  which  is  defined  on  a  positive  range  and 
bounded  will  be  used  to  model  the  initial  enzyme  content. 

The  beta  probability  density  family  is  represented  by  a  two-parameter  density 
function  typically  defined  over  the  closed  interval  between  zero  and  one  [12]  and  is  given 
by 


/w= 


Tja  +  P) 

0, 


X 


(l-x) 


p-\ 


0  <  X  <  1,  a  >0,  >  0, 

X  >  1. 


Here  r(a)  is  the  gamma  function  [12]  where 


r(cjr)  =  x“  ""dx. 

Refer  to  Ref.  12  for  additional  properties  of  the  gamma  function. 
The  beta  cumulative  distribution  is  then 


(4.9) 
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P(x<Xf^)  = 


T(a  +  J3) 


r(cc)r(j3) 


'(1-x)'^  ^dx. 


0  <  Xq  <  1 


(4.10) 


[1,  Xq  >  1 

This  beta  distribution  can  be  expressed  in  terms  of  the  incomplete  beta  function 
given  by 


Thus, 


(4.11) 


P(x<Xo)  = 


T(a  +  J3) 

r(a)n/3) 


0  <  Xq  <  1 


[1,  Xq  >  1 

As  noted  above,  the  concentration  of  each  enzyme  within  the  spore  is  between 
zero  and  where£^(r,4)  =  ^^/^(r,^) .  £®^^is  the  maximum 

concentration  of  each  enzyme  prior  to  any  treatment  and/(r,r^)  is  given  in  Equation  4.5. 
By  dividing  the  concentration  of  each  enzyme  amount  by£'^‘^(r,t^),  the  concentration  of 

each  enzyme  will  be  proportions  between  zero  and  one.  In  order  to  generate  the 
physically  realistic  representation  of  the  spore’s  initial  enzyme  content,  it  was  assumed 
that  initially  a»  P .  From  this  assumption,  the  mean  amount  of  enzymes  will  be  close  to 
the  maximum  possible,  Emax,  and  the  probability  distribution  of  enzymes  will  be  skewed. 
Although,  the  beta  density  function  is  only  nonzero  on  the  interval  0  <  x  <  1,  it 


can  be  rescaled  for  an  arbitrary  finite  interval.  To  this  end,  setx  = 


a 


where 


a  =  E^(T,tp  and  define  E^  =  XqA  in  Equation  4.10  to  produce 
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This  leads  to  the  associated  scaled  density  function 

£•  {t  )  (a  -  £  (t  )y-' ,  0<E(n<a 

nE^(0)=\na)n/^)  (4.12) 

0,  E^{t^)>a 

Figure  4.3  shows  an  illustration  of  a  notional  enzyme  fitness  density  function  with 
a  =  10  and  /?  =  3  prior  to  any  heat  treatment.  The  random  variable,  Eq,  will  be  drawn  from 


this  enzyme  fitness  distribution  with  mean,  and  variance,  <7^^ .  As  exposure  time 


increases,  the  distribution  of  the  enzymes  will  evolve.  Since  Eq  is  a  random  variable, 
Er(te)  is  a  random  variable  with  mean,  and  variance,  ^ . 
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Figure  4.3:  Enzyme  Fitness  Prior  to  Any  Heating  (Beta) 

Using  the  expectation  rules  above, 


=  E{EJ,{T,0} 
=  f,{T,OE[E,} 
~  fs  K  )/^£o 
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=  E{E,\f,{T,Qf]- 
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Using,  the  property  that  xr(x)  =  r(x  +  l) 

Yi^cc -\- P -\-Y)  cc  (cr  +  l)r(cr + 1) 


cTE:=vfE(T,oY 


=  UEiT,Of 


=  UEiT,Of 


^  a 


{cc^ Y(^oc -yV)  i^oc -\- P -\-V)Y(^(x -\- P -\-V)  ^oc-yp 
a{a  +  Y) 


{a  +  P){a  +  P  +  Y) 

aP 


^  a 
oc^  P 


{a  +  pf{a  +  P  +  \) 
=  UE(r,Ofal 


As  exposure  time  inereases  from  =  0  to  =T2,  Er(te)  moves  eloser  to  the  kill 


threshold  as  depicted  in  Figure  4.4.  The  curve  outlined  in  blue  represents  the  enzyme 
distribution  at  exposure  time,  =  0.  Since  no  enzyme  damage  has  yet  occurred,  the 


population  mean  is  10/13  (for  a=10  and  |3=3).  The  curve  outlined  in  green  represents  the 
enzyme  distribution  after  some  enzyme  damage  at  time,  4  =  Tj  ,  Tj  <  Tj  .  Now  that  damage 
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has  been  done  to  the  enzymes,  the  mean//^.  (Tj)  < //^^ .  The  curve  outlined  in  red 
represents  the  enzyme  distribution  at  exposure  time, =  Tj. Here  the  spore’s  enzymes 
have  accumulated  even  greater  damage  and  the  mean  jue  (^2)  <  Also,  the  majority  of 

the  fitness  distribution  has  exceeded  the  kill  threshold  and  thus,  there  will  be  little  viable 
enzymes  available  for  DNA  repair  process.  Possessing  behaviors  similar  to  those  when 
using  the  Gaussian  representation,  as  exposure  time  increases  the  mean  and  variance  of 
the  distribution  also  decreases. 


Figure  4.4:  Enzyme  Fitness  Evolution  (Beta) 

The  probability  of  kill,  Pk,  can  be  found  by  computing  the  cumulative  probability 
up  to  the  kill  threshold.  Notationally, 
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< 


yE^(T,t^)  E^{T,Qj 


=  P,  EJ(T,tJ<E, 


k  ’ 


~  E/c  Eq  <  E^. 


where  E, 


f(T,0 


Then  from  equation  4.10, 


P,iE,<E,)  = 


Tia  +  ^)  A 


Tia)r(^) 


l'E,'^-\l-E,y-^dE„  0<E,<1 


1, 


E,>1 


r(a  +  /]) 

1, 


B,  {a,j3). 


0  <  £■,  <  1 


E,>1 


Probability  of  Survival. 
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Figure  4.5:  Gaussian  and  Beta  Distribution  Comparison 
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Figure  4.5  shows  a  comparison  between  the  Gaussian  and  beta  density  functions 


prior  to  any  exposure  time.  In  order  to  make  this  comparison,  the  mean  and  variance  were 
chosen  to  be  the  same,  that  is, 

~  0.7692, 
o-/  «  0.0127. 

The  Gaussian  density  function  exceeds  one  while  the  beta  density  function  is  bounded  by 
one.  This  will  be  the  basis  for  comparing  the  probability  of  survival  for  the  Gaussian  and 
beta  distribution. 

The  notional  parameters  used  within  the  probability  of  kill  models  can  be  found  in 
Table  4.1  [5,22,  25,29,  40]. 


Table  4.1:  Parameters  for  Probability  of  Kill  Models 


Parameter 

Value 

Units 

R 

1.986  X  10'^ 

kcal/morK 

A; 

4x  10'" 

sec'' 

El 

26.1 

kcal/mol 

A2 

4.828  X  10'= 

sec'' 

E2 

38.3 

kcal/mol 

A 

1.6  X  10^ 

sec'' 

E, 

10 

kcal/mol 

Figure  4.6  depicts  the  probability  of  enzyme  survival  during  dry  heating  or7  =  0 
with  an  initial  release  temperature  of  225°C  and  a  kill  threshold  of  0.1.  That  is,  once  the 
enzyme  population  drops  to  10%,  the  spore  will  not  contain  the  necessary  enzyme  for 
DNA  repair.  The  red  curve  is  the  beta  distribution  representation  for  the  fitness  model 
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and  blue  curve  is  the  Gaussian  representation.  As  it  can  be  seen,  the  beta  representation 
drops  off  faster  than  the  Gaussian.  This  is  because  the  Gaussian  distribution  is  not 
bounded  above  like  the  beta  distribution.  With  the  beta  representation,  it  will  take  about 
thirty  seconds  to  achieve  a  significant  reduction  (lO'^*^)  in  enzyme  survival  and  for  the 
Gaussian  representation,  it  will  take  about  thirty-six  seconds  to  receive  the  same 
reduction. 


Figure  4.6:  Probability  of  Enzyme  Survival  (225  C) 

Figure  4.7  illustrates  the  probability  of  enzyme  survival  for  dry  heating  with  an 
initial  release  temperature  of  300°C.  Again,  notice  that  the  estimate  of  spore  survival  in 
the  beta  model  falls  of  faster  than  the  Gaussian  model.  When  heated  at  higher 
temperatures,  equivalent  enzyme  damage  takes  place  within  a  shorter  duration.  With  the 
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beta  model,  it  only  takes  approximately  0.0038  seconds  to  reduce  the  enzyme  survival  by 
a  significant  amount.  With  the  Gaussian  model,  it  takes  approximately  0.0041  seconds  to 
receive  the  same  reduction.  This  can  be  explained  by  the  fact  that  from  Chapter  II  at 
290°C,  chemically  bound  water  becomes  available  within  the  spores.  The  release  of 
chemically  bound  water  causes  the  enzymes  to  be  damaged  faster  as  more  water  is 
available  for  hydrolysis  reactions. 


Figure  4.7:  Probability  of  Enzyme  Survival  (300  C) 
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V.  Germination  and  Outgrowth 


Introduction 

Once  the  environment  becomes  nutrient  rich,  the  spore  can  break  the  dormant 
state  and  move  into  germination.  There  are  three  stages  of  germination:  activation, 
germination,  and  outgrowth.  During  activation,  the  spore  prepares  to  germinate.  During 
germination,  the  dormant  spore’s  resistance  mechanisms  are  eliminated  and  metabolic 
activity  begins  as  a  vegetative  cell  is  formed.  Also,  during  germination,  the  spore’s  DNA 
will  be  repaired  if  there  are  enough  viable  enzymes  for  the  repair  process.  If  the  spore’s 
DNA  is  repaired,  the  spore  transitions  into  outgrowth.  During  outgrowth,  the  spore’s 
proteins  and  DNA  are  synthesized  for  daughter  cells.  An  initial  vegetative  cell  is  first 
formed  from  the  spore.  Whether  the  spore  can  survive  to  the  last  stage  of  outgrowth  is 
dependent  on  several  events,  such  as  the  exposure  time  and  temperature  during  thermal 
inactivation,  the  initial  state  of  the  DNA,  as  well  as  the  extent  of  viable  enzymes  to 
complete  the  repair  process. 

Germination 

Germination  is  an  irreversible  and  degradative  process  in  which  the  spore 
transforms  from  dormant  and  highly  resistance  spore  to  an  actively  growing  and  dividing 
cell.  Many  changes  occur  during  germination.  The  spore’  resistance  to  heat,  desiccation, 
ultraviolet  radiation,  and  chemicals  are  lost. 

Activation  reversibly  conditions  the  spore  to  germinate  if  appropriate  germinants 
are  provided  [27].  This  process  may  take  less  than  one  minute  for  an  individual  spore  and 


50 


longer  for  a  population  of  spores.  Mild  heating  of  approximately  65 °C  or  the  presence  of 
certain  chemicals  triggers  this  stage.  After  activation,  “the  core  enzymes  do  not  function 
properly  and  the  core  water  content  remains  low”  [27]. 

Germination  terminates  the  dormant  spore  resistance,  metabolic  inactivity,  and 
refractivity.  It  occurs  when  the  spore  come  into  contact  with  many  types  of  nutrients  such 
as  amino  acids,  sugars,  and  purine  nucleosides  [33].  The  result  is  the  initiation  of  a 
breakdown  of  macromolecules  and  excretion  of  spore  substances  such  as  the  release  of 
the  dipicolinic  acid,  calcium,  and  peptides  [16].  These  substances  are  replaced  by  water 
[30].  There  are  observable  physical  changes  within  the  layers  of  the  spore  structure 
during  this  stage.  The  spore  coat  remains  intact,  but  the  cortex  dissolves  as  the  core 
swells,  filling  the  exosporium.  This  resulting  increase  in  the  spore  volume  is  an  important 
step  in  the  rehydration  of  the  core. 

Furthermore,  during  this  stage,  the  spore’s  damaged  DNA  is  repaired.  As 
mentioned  above,  since  enzyme  synthesis  does  not  occur  during  dormancy,  the  spore 
must  contain  enough  viable  enzymes  in  order  to  complete  the  DNA  repair  process. 
Additionally,  “if  too  much  damage  has  been  accumulated  during  spore  dormancy,  this 
damage  can  overwhelm  the  capacity  of  the  repair  systems  and  lead  to  the  death  of  the 
germinated  spore”  [34]. 

As  mentioned  in  Chapter  IV,  although  there  are  numerous  types  of  mechanisms  to 
repair  damaged  DNA,  it  is  believed  that  BER  will  mostly  likely  be  necessary  within  the 
spore  during  germination. 
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DNA  BER. 


Figure  5.1  illustrates  the  BER  process  after  U  has  been  formed  by  deamination  of 
C  and  is  therefore  opposite  G  in  the  DNA  strand.  Initially,  the  bond  between  uracil  and 
the  deoxyribose  is  cleaved  by  a  DNA  glycosylase.  This  leaves  an  apurinic/apyrimidinic 
(AP)  site,  a  sugar  with  no  base  attached  in  the  DNA.  Next,  AP  endonuclease  cleaves  the 
DNA  chain.  Deoxyribose  phosphodiesterase  removes  the  remaining  deoxyribose.  The 
gap  is  then  filled  by  DNA  polymerase  and  sealed  by  ligase  with  the  correct  base, 
cytosine. 


DNA  contuininy  U  formed  by 
de-Jiniruition  ol  C 


Al’sile 


/  ''  \ 

J 

n 

1 

B 

i 

vV 

\ 

Figure  5.1:  Base  Excision  Repair  [3] 
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Many  pathways  exist  for  the  BER.  Quantitative  models  exist  for  this  type  of 
repair  mechanism  [36].  These  models  were  based  on  Michaelis-Menten  kinetics,  which  is 
also  the  foundation  of  the  work  in  this  research 
Michaelis-Menten  Kinetics 


Enzymes  act  as  regulators  of  biochemical  processes  such  as  DNA  repair.  Enzyme 
kinetics  is  the  study  of  the  rates  of  enzyme-catalyzed  reactions.  Michaelis-Menten 
kinetics  is  a  model  of  enzyme  kinetics.  The  Michaelis-Menten  equation  describes  the 
rates  of  enzymatic  reactions  by  relating  reaction  rate  to  the  substrate  concentration.  The 
main  assumption  is  that  the  enzyme  and  the  substrate  react  reversibly  to  initially  form  a 
complex.  The  complex  then  breaks  down  to  form  the  free  enzyme  plus  one  or  more 
products,  as  shown  below. 


Equation  5.1  states  that  an  enzyme  E  combines  with  substrate  S  to  form  an  ES  complex, 
with  a  rate  constant  of  ki.  The  ES  complex  then  can  dissociate  to  E  and  S,  with  a  rate 
constant  k.i  or  it  can  form  a  product  P,  with  a  rate  constant  ka.  The  turnover  number,  kcat, 
is  the  number  of  substrate  molecules  converted  into  product  by  an  enzyme  molecule 
when  the  enzyme  is  fully  saturated  with  substrate  and  is  equal  to  the  rate  constant,  k2. 

Using  Equation  5.1,  rate  equations  can  be  formed  for  the  concentration  of 
enzyme,  substrate,  complex,  and  product,  (the  notation  [  ]  represents  concentration  of) 

^  =  UES^,  (5.2) 

dt 


J[5] 

dt 


(5.3) 
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=  k^AE][S]-(k,  +k_,)[ES] , 


(5.4) 


^  =  -kJE]  [5]  +  (k,  +  k_,  )[£5] ,  (5.5) 

dt 

where  t  represents  time. 

The  velocity  of  the  reaction,  V,  is  the  rate  of  product  of  the  concentration  of  a 
substrate  with  the  rate  coefficient.  For  example  in  Equation  5.2,  the  reaction  velocity  for 
[P]  is  the  product  of  rate  constant,  ka,  with  the  concentration  of  complex,  [ES],  i.e. 

V  =  ^^  =  UES].  (5.6) 

dt^ 

To  solve  Equations  5. 2-5.6,  initial  conditions  are  necessary.  If  only  the  enzyme  and 
substrate  are  present  at  time  t=0,  then 


([SI  [E],  [E5],  [P])  =  (5o,  E„  0, 0)  (5.7) 

Briggs  and  Haldane  suggested  the  steady  state  for  the  complex  is  when  the  concentration 
of  complex  stays  the  same  even  if  the  concentrations  of  starting  materials  and  products 
are  changing.  That  is,  the  rate  of  change  in  the  complex,  [ES],  is  equal  to  zero,  i.e. 

k^,  [P]  [S]  -  (k^  +  k_,  )[P5]  =  0. 


Therefore, 


[E][S]_(k,+k_,) 


Equation  5.8  is  the  Michaelis  constant.  Km, 


(5.8) 


Km  = 


(^2+^-l) 


(5.9) 


The  Michaelis  constant  is  independent  of  enzyme  and  substrate  concentrations.  Also,  by 
adding  Equations  5.4  and  5.5, 
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This  implies  that 


4  [ES]  +  [E]  =0. 
dt 


m+[ES}=c, 

where  C  is  a  constant.  From  Equation  5.7,  at  t=0, 


[E]  +  [E5]  =  Eo 

Now,  by  substituting  Equation  5.10  into  Equation  5.8, 


and  solving  for  [ES], 


(5.10) 


[E5]=  . 


Therefore,  from  Equation  5.6, 


k2lS]EQ 


(5.11) 


(5.12) 


This  gives  the  Michaelis-Menten  equation.  At  low  substrate  concentration,  when  [S]  is 
much  less  than  Km, 


k2lS]EQ 

Km 


(5.13) 


The  reaction  becomes  first  order  with  the  rate  directly  proportional  to  the  substrate 
concentration.  At  high  substrate  concentration,  when  [S]  is  much  greater  than  Km, 


y  =  ^2^o-  (5-14) 

The  reaction  is  now  independent  of  substrate  concentration  and  the  rate  becomes 
maximal.  The  reaction  velocity  is  then  said  to  saturate  with  respect  to  the  substrate  or  to 
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exhibit  saturation.  Using  the  Michaelis-Menten  theory,  quantitative  models  were 


developed  for  the  BER  of  DNA  [36].  These  models  will  be  further  explored  in  the  next 
section. 


DNA  BER  Model 


Figure  5.2  shows  an  illustration  of  a  pathway  for  the  BER  mechanism  with  the 
definitions  of  substrate  concentrations  yi  (nM),  enzyme  concentrations  Ci  (nM),  and 
reaction  velocities  Vi  (nM/sec). 


y2 


yi 


Abasic  Site 


Gapped, 

5dRp 


ya 

Nicked, 

5dRp 


ys 


Nicked 


Repaired 
DNA  Adduct 


Apei 


PolB, 
gap  filling 


Pol|i, 

dRP  Lyase 


LigI 


Figure  5.2:  Schematic  of  the  DNA  Repair  Model  (BER)  [36] 

Initially,  Apel  incises  the  abasic  site  producing  a  single  nucleotide  gap  with  a  3’- 
hydroxyl  group  and  a  5’-dRp  group.  Then,  DNA  Poip  (e2)  fills  the  gap  producing  a 
nicked  substrate  with  a  5’-dRp  remaining.  The  5’-terminus  is  then  processed  by  the 
dRp  lyase  activity  of  Poip.  Lastly,  DNA  ligase  (e3)  seals  the  nicked  substrate. 


A  quantitative  model  of  human  DNA  BER  was  developed  from  kinetic  data 
available  in  by  Sokhansanj  et  al.  [36].  By  assuming  Michaelis-Menten  kinetics,  the  model 
equations  are 
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dt^ 


=  -Vi, 


dt. 


=  v,  -v,, 


dt. 


=  v,-v.. 


dt. 


=  v,-v.. 


Vi  = 


^  kA  ^ 

3^1  +  ^i 


^2  =>’2 


V3  =>’3 


V4=>’4 


^  ^2^2  ^ 

>’2  +  ^2  y 
A 


3  y 


^4^3  ^ 


4  J 


where  k  (1/sec)  represents  kcat,  K  (nM)  represents  the  Michaelis  constant,  and  h 
represents  the  repair  time  (sec).  Notice  that  the  velocity  rates  are  dependent  on  the 
individual  enzyme  concentrations.  Since  the  enzymes  too  are  damaged  during  the 
inactivation  process,  the  amount  of  viable  enzymes  is  decreased.  If  the  spore  contains  a 
sufficient  amount  of  enzymes,  then  the  repair  time  will  be  quick.  Whereas,  if  the  spore’s 
enzymes  are  reduced  substantially,  the  repair  time  will  take  longer. 

Each  compartment  in  Figure  5.2  is  dependent  on  the  previous  one.  If  ei  is  reduced 
during  the  thermal  inactivation,  then  vi  would  be  slow  and  it  will  affect  the  speed  of  the 
second  compartment,  y2.  This  in  turn  would  affect  the  third  compartment,  ys.  And  so  on. 
The  same  can  be  same  for  02  and  03.  Reductions  in  the  enzyme  concentrations  will  also 
affect  the  repair  velocities  of  each  compartment. 

This  model  has  been  mimicked  to  explain  the  repair  process  within  the  Bacillus 
spores.  Furthermore,  this  research  is  interested  in  DNA  information  content  instead  of 
DNA  concentration.  Therefore,  a  dimensionless  variable  is  introduced.  Also,  from 
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Chapter  IV  since  the  fitness  of  each  enzyme  is  assumed  to  come  from  a  beta  distribution, 
the  enzyme  concentrations  will  be  scaled  by  .  This  results  in  the  following 


e,iT,0  =  r,EZ(T,0 


(5.15) 


where  represents  the  maximum  enzyme  concentration  for  the  i*  enzyme  type 

after  exposure  time.  Further,  Iq  represents  the  initial  state  of  the  total  DNA  (critical  and 
non  critical)  such  that 


Ic=hO^ 


and  Ic  represents  the  initial  state  of  critical  DNA  with  0  the  fraction  of  the  total  DNA 
which  is  critical.  Moreover,  the  repair  process  does  not  recognize  the  difference  between 
critical  and  non  critical  DNA  and  therefore  will  repair  the  total  damaged  DNA. 
Substituting  the  dimensionless  variables  into  the  DNA  BER  model  produces  five 
equations  as  follows 


dt^ 

d^2(0 

dt^ 

d<^2(tr) 

dt^ 

d^AO 

dt^ 

d^AK) 

dt^ 


=  -^Atr) 
=  ^Atr) 

=  ^2iO 

=  ^Atr) 
=  ^Atr) 


k,r,El(T.t,) 


KnEZiT.Q 

V  4fi,)h+K,  J 

[  4AK)h+K,  J 

k,r2E‘S,(T,t,) 


&(0 


KnEZiT.t.) 

I  / 

KnE^ZiT.O 

(  ^,((,)/o+a:3  J' 

KriEZiTA) 


^AOi,+K. 


4  J 


KnEZ(T,0 
[  ^AOh  +  K,  , 


(5.16) 


In  this  model,  when,^.(t^)/Q  »K.,  the  rates  become  independent  of  the  substrate 


concentration  as  in  Equation  5.14.  When^.(t^)/^  «  K.,  the  rates  become  first  order  linear 
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equations  as  in  Equation  5.13.  The  initial  concentration  of  DNA,  Iq,  plays  an  important 
role  in  the  rate  equations.  To  further  understand  the  influences  of  Iq  on  the  DNA  repair 
process,  a  sensitivity  analysis  was  conducted. 

Sensitivity  Analysis. 

Table  5.1  contains  the  values  for  parameters  used  within  the  original  DNA  BER 
model  [36]. 


Table  5.1:  Parameters  for  Model 


Parameter 

Value 

Parameter 

Value 

ki 

0.625  sec  ' 

Ki 

300  nM 

k2 

0.04  sec  ' 

K2 

100  nM 

ks 

0.78  sec'' 

K3 

39  nM 

U 

0.05  sec  ' 

K4 

130  nM 

Additionally,  from  Chapter  IV  the  DNA  information  content  and  enzymes 
diminishes  exponentially  from  hydrolysis  and  pyrolysis  damage,  i.e.,  recall 

Io(0  =  IJd(T,0, 
e,(0  =  ej,(t,o. 

If  the  rate  coefficients  associated  with  pyrolysis  and  hydrolysis  for  both  DNA 
information  content  and  enzymes  are  the  same,  then  the  spore’s  DNA  and  enzymes  will 
be  reduced  at  the  same  rate.  For  example  by  this  simplification,  after  the  inactivation 
process,  if  10%  of  the  spore’s  DNA  is  damaged  then  10%  of  the  spore’s  enzymes  are  also 
damaged. 
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An  additional  assumption  must  be  made  about  the  maximum  enzyme 
concentration  for  each  enzyme  prior  to  thermal  heating.  The  spore  contains  an  initial 
concentration  of  enzymes  used  for  the  repair  process.  These  initial  maximum 
concentrations  are  presented  in  Table  5.2.  These  values  were  chosen  based  on  the 
parameter  values  in  Ref.  33,  an  exposure  time  of  2  seconds,  and  temperature  of  300°C. 

Table  5.2:  Enzyme  Maximum  Parameters 


Enzyme 

Rd) 

maxQ 

e'Kt.i.) 

l(Apel) 

130  nM 

83.19  nM 

2  (Poip) 

50  nM 

32.00  nM 

3  (Ligl) 

150  nM 

95.99  nM 

Multiple  cases  were  examined.  Case  1  assumed  «K.  for  all/  =  l,...,4. 

Case  2  assumed  x  for  some/ =  1,..., 4.  Case  3  assumed  K.  for 

all/  =  l,...,4.  The  initial  DNA  information  parameter,  lo,  for  each  of  these  cases  is  given 
in  Table  5.3, 


Table  5.3:  lo  values  for  Cases  1, 2, 3 


Case 

lo 

1 

10  nM 

2 

100  nM 

3 

1000  nM 
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Further,  three  levels  of  damaged  DNA  and  enzymes  were  considered.  These  were 
when  10%,  50%,  and  90%  of  DNA  and  enzymes  are  damaged  during  the  thermal  heating. 
Figure  5.3  illustrates  the  repair  of  DNA  when  10%  of  the  spore’s  DNA  and  enzymes  are 
damaged.  In  this  instance,  there  is  little  difference  between  lo  =  10  nM  and  lo  =  100  nM. 
However,  for  Iq  =  1000  nM,  the  rate  of  repair  decreases  and  it  takes  longer  to  repair  the 
damaged  DNA.  For  example,  it  takes  approximately  100  seconds  to  repair  half  of  the 
damaged  DNA  whenl^  =  10  nM  and  approximately  146  seconds  when  lo  =  1000  nM. 


Figure  5.3:  DNA  Repair  Analysis  when  10%  of  DNA  &  Enzymes  are  Damaged 

Figure  5.4  depicts  the  repair  of  DNA  when  50%  of  the  spore’s  DNA  and  enzymes 
have  been  damaged.  Compared  to  Figure  5.3,  the  repair  time  has  increased  because  lower 
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enzymes  concentrations  are  now  available  for  the  repair  process.  Again,  when  Iq  =  10  nM 
and  lo  =  100  nM,  there  is  not  great  difference  in  the  time  required  to  repair  the  damaged 
DNA.  It  takes  about  187  seconds  whenlg  =10  nM  and  about  223  seconds  when  Iq  =  100 

nM  to  repair  half  of  the  damaged  DNA.  Whereas,  for  Iq  =  1000  nM,  it  takes 
approximately  580  seconds  to  repair  half  of  the  damaged  DNA. 


Figure  5.4:  DNA  Repair  Analysis  when  50%  of  DNA  &  Enzymes  are  Damaged 

Figure  5.5  illustrates  the  repair  of  DNA  when  90%  of  the  spore’s  DNA  and 
enzymes  have  been  damaged.  When  the  majority  of  the  spore’s  enzymes  are  damaged, 
the  repair  process  takes  much  longer  than  in  the  previous  cases.  In  order  to  repair  half  of 


62 


the  damaged  DNA  when  Iq  =  1000  nM  ,  the  time  required  would  be  4477  seconds  or 


roughly  1  hour  and  14  minutes. 


Figure  5.5:  DNA  Repair  Analysis  when  90%  of  DNA  &  Enzymes  are  Damaged 

Thus,  the  initial  state  of  DNA  is  important  in  the  repair  rate  and  whether  or  not 
the  spore  will  be  able  to  repair  the  damaged  DNA.  If  Case  1  holds  true,  then  the  DNA 
BER  model  becomes  a  system  of  first  order  linear  equations.  Therefore,  an  analytical 
solution  can  be  written.  The  influence  that  assumption  has  on  the  DNA  BER  model  will 
be  examined  in  the  next  section. 
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Linear  Model  ofDNA  BER  Model. 

When^aj/,  «  ,  the  substrates  stay  far  from  saturation  where  a  saturation 


value  will  be  defined  as  .  The  DNA  BER  model  equations  can  be  reduced  to  a  first 

order  linear  system.  Assume  that  «K^  for  i=l,...,4.  The  influence  of  the 

enzymes  and  initial  state  of  DNA  will  be  explored. 

To  linearize  the  nonlinear  model  in  Equation  5.16,  the  Taylor  series  expansion  of 
functions  was  used.  Let/f^j  be  a  function  such  that/(0i=0  and  define  the  differential 
^iitr)  =  fidit^)) .  If  ^(d)  =  0,  fi(^)  may  be  expanded  in  a  Taylor  series  about  zero  to 
produce 

/;.(|)  =  /;.(|(0))+y/;.(|(0)).(|-|(0))+o((|-|(0))"),  / =i,...,4. 


Then  since /;.(^(0))  =  /;.(0)  =  0, 


Using  Equation  5.16,  the  following  is  obtained 


/  =  !,...,  4. 


v/.(l)U=- 


V/2(|)Lo=- 


(1) 


r(i) 


-,0,0,0  , 


(2) 


v/3(l)U=  0, 


K,  K, 

klYl^mca  k^Yl^max 


-,o,o  , 


K,  K, 

7(2) 


,0 


r(3) 

'^max 


Vf(c5)l  =Jo  0  _KY3^ 

In  these  results  the  exposure  time  and  temperature  dependence  on£'^'^(r,r^)  has  been 


suppressed. 
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For  simplicity,  let 


.  KnEZc  2  3  2 

A-  ^  ^  ^  ’A;- 

XVj^  “^^3  ^^4 


(5.17) 


Then  the  Jacobian  for  the  system  beeomes 


Thus,  =  J^{t^).  Therefore,  a  simple  solution  can  be  given  as  =  e'’*' Since  J  is 
a  constant  matrix,  {-/IJti  distinct  eigenvalues  of  J  with  corresponding 


eigen veetors  l^^ithe  solution  can  be  given  as 


m=i=Yc,i,=z. 


where  Z  is  the  matrix  of  eigenvectors. 


.c=z-‘4 


Z  =  [4'„C2,C3,4'4]- 


The  eigenvalues  of  J  are-Zj,  -Z4  with  corresponding  eigenvectors. 


Ci=  1,- 


\  X2 


A^—  X2  A^—  Au2  A^—  A^  A^—  A^  A^—  A^  A^—  Az^ 


^2  =  0,  1,  - 


^  ^ _ ^ 

A^~  A[^  A^~  A^  A^~  Ai\ 
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^3  = 


0,  0,  1, 


-[T 


C4  =  0,  0,  0,  1  . 


Figure  5.6  shows  comparisons  between  the  numerical  solution  of  the  nonlinear 
DNA  BER  model  equations  and  the  linearization  solution  if  50%  of  the  spore  DNA  and 
enzymes  are  damaged  and  Io=10  nM.  Values  from  Tables  5. 1-5.2  were  used  with  an 
exposure  time  of  2  seconds  at  a  temperature  of  SOO^C.  For  and  ^2,  there  is  no 
distinction  between  the  linear  and  nonlinear  equations.  For  ^3  and  ^4  ,  the  linear  and 
nonlinear  equations  reach  different  peaks  and  as  time  progresses  they  evolve  together. 
Given  the  Michaelis  constant,  Ki,  in  Table  5.1  and  lo  =  10  nM,  the  saturation  value  for  ^1, 
^2,  ^3,  ^4  respectively  are  30,  10,  3.9,  and  13.  Shown  in  Figure  5.6,  each  stays  far  from 
its  respective  saturation  value. 


Figure  5.6:  Nonlinear  and  Linear  Model  Comparison  (Io=10  nM) 
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Figure  5.7  shows  comparisons  between  the  nonlinear  DNA  BER  model  and  the 
linear  model  when  Io=100  nM  and  50%  of  the  spore  DNA  and  enzymes  are  damaged.  For 
^1,  there  is  again  no  distinction  between  the  linear  and  nonlinear  model  equations.  For  ^2, 
as  exposure  time  progress  there  is  some  difference  between  the  nonlinear  and  linear 
equations.  For  ^3  and  ^4 ,  the  linear  and  nonlinear  equations  reach  different  peaks  but  this 
difference  is  not  significant  since  the  scaling  is  small.  Further,  as  repair  time  progresses, 
the  linear  and  nonlinear  equations  evolve  together.  Here  the  saturation  values  for  ^2,  ^3, 
^4  respectively  are  3,  1,  0.39,  and  1.3.  Each  stays  away  from  its  respective  saturation 
value. 


Figure  5.7:  Nonlinear  and  Linear  Model  Comparison  (Io=100  nM) 
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Figure  5.8  shows  comparisons  between  the  nonlinear  DNA  BER  model  and  the 
linear  model  when  Io=1000  nM  and  50%  of  the  spore  DNA  and  enzymes  are  damaged. 
For  ^1,  early  on  there  is  a  difference  between  the  linear  and  nonlinear  model  equations  but 
as  repair  time  progress,  they  evolve  together.  For  ^2,  the  linear  and  nonlinear  equations 
reach  different  peaks  and  as  exposure  time  progress  there  is  some  difference  between  the 
two.  For  ^3  and  ^4 ,  the  linear  and  nonlinear  equations  reach  different  peaks.  There  is  also 
a  point  at  which  the  nonlinear  equations  remain  constant  and  therefore,  the  nonlinear  and 
linear  equations  do  not  evolve  together.  When  Io=1000  nM,  the  saturation  values  for  ^1, 
^2,  ^3,  ^4  respectively  are  0.3,  0.1,  0.039,  and  0.13.  For  and  ^2,  each  passes  their 
respective  saturation  value.  For  ^3,  the  saturation  value  is  not  reach  and  for  ^4,  the 
saturation  value  is  approached. 


Figure  5.8:  Nonlinear  and  Linear  Model  Comparison  (Io=1000  nM) 
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The  repair  term  in  the  DNA  repair  model  can  be  given  as 


=  I'  4(r) 


''  Kr,EL  ' 

dr. 


dr. 


K,+Ur)h 


oj 


Assume  «i<r4,  then 


6^.)  = 


KYiK 

K. 


£  ^4(.T^)dT. 


Evaluating  the  integral  in  Equation  5.18  produces 


^5(^r;r)  =  ^r^^ei 


+—  \-e 


+ 


(5.18) 


(5.19) 


where 


e,  = 


\ _ ^ _ ^ 

x^— ^2  X^— X^  X^— Xj^ 

X^  X2  x^ 
A^—  X2  X2~ X^  X2  ~ X^ 

A 

X^—  X2  x^  —  x^  x^  —  x^ 


Observe  that  under  the  assumption  each  enzyme  is  subjected  to  the  same  rate  of 
hydrolysis  and  pyrolysis  reactions, /^(r, 4)  can  be  factored  with  respect  to  resulting  in 
the  exposure  time  appearing  only  in  the  exponential  terms  in  Equation  5.19.  Therefore, 
the  effect  of  exposure  time  will  only  change  the  rate  at  which  the  velocity  equilibrates. 

If  the  spore’s  DNA  is  repaired  without  complications,  the  next  step  in  germination 
is  outgrowth. 
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Outgrowth 

Outgrowth  is  the  development  of  a  vegetative  cell  from  a  germinated  spore  [17]. 
Outgrowth  is  dependent  on  the  repair  of  the  protein-synthesizing  system  and  an  ordered 
synthesis  of  proteins  [21],  According  to  Hansen  et  al,  the  synthesis  of  RNA  begins  after 
germination  and  then  the  proteins  are  synthesized  when  the  new  vegetative  cell  emerges 
and  divides  and  outgrowth  can  last  for  hours  in  one  spore,  but  in  a  population  of  spores,  it 
occurs  synchronously  [17].  Let  outgrowth  time,  to,  be  defined  as  the  time  when  the 
experimental  data  is  taken  to  determine  if  the  spore  produced  a  viable  cell. 

During  outgrowth,  the  core  wall  thickens  to  become  the  vegetative  cell  wall. 
Remains  of  the  outer  part  of  the  old  spore  coat  and  the  exosporium  are  found  within  the 
new  vegetative  cell  [16].  Also,  in  the  vegetative  cell,  “the  heat  resistance  is  increased  and 
hardiness  is  gone”  [39].  Figure  5.9  illustrates  the  vegetative  cells  forming  from  the 
spores. 

The  spore’s  chance  of  survival  is  dependent  on  numerous  environmental 
parameters.  The  probability  of  survival  models  examined  in  Chapter  IV  was  based  on 
knowledge  of  a  critical  threshold  for  the  enzymes.  In  Figures  4.6  and  4.7,  this  threshold 
was  chosen  to  be  0.1.  That  is,  once  the  enzyme  population  drops  to  10%,  the  spore  will 
not  contain  the  necessary  enzyme  for  DNA  repair.  However,  this  threshold  is  not  a 
known  quantity.  In  fact,  if  no  DNA  damage  occurs,  no  repair  enzymes  are  necessary  for 
outgrowth.  Thus,  the  enzyme  threshold  may  also  depend  on  DNA  damage.  The  influence 
that  the  DNA  threshold,  outgrowth  time,  and  viable  enzymes  have  on  the  spore’s 
probability  of  survival  will  be  examined  in  the  following  section. 
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Figure  5.9:  Illustration  of  Bacillus  subtilis  Vegetative  Cells  [9] 

The  image  is  of  a  stained  preparation  of  Bacillus  subtilis  showing  endospores  as 

green  and  the  vegetative  cell  as  red. 


Probability  of  Survival 

Multivariate  Probability  Distribution. 

Let  Yi,  Y2,  Y„  be  continuous  random  variables  with  joint  cumulative 
distribution  function  F  (yj,  y2,  y„).  That  is, 

F{y,,y^,...,yn)=  T  P  •••  f{t,,f,...,tjdf...di,  (5.20) 

J—co  J  —00  J—co 

where  the  function/(jj,  the  joint  probability  density  function  [12]. 

If  Fi,  F2, ...,  Yn  are  independent  then 


F{y„y2,...,y„)  =  F,  {y,  (y^)  •  •  •F„  (y„ ) 


and  the  joint  density  function  becomes 
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/(yi,y2,-,y„)  =  /i(yi)/2(y2)"-f„(y„)[i2]. 

Using  these  concepts,  the  probability  of  survival  model  based  only  on  single  DNA 
threshold  will  be  derived. 

Enzyme  Dependence. 

The  spore’s  chance  of  survival  is  dependent  on  the  available  viable  enzymes  for 
the  DNA  repair  process.  Therefore,  the  total  number  of  possible  combinations  of  each 
enzyme,  Ci,  ca,  and  es  will  determine  the  domain  in  which  the  spore  has  a  chance  of 
survival.  Assuming  these  enzymes  are  independent  random  variables,  the  probability  of 
kill  is  defined  as 


P, 


yj  t^,e  +i^{T,te)<yc 


=  III  n  /,(«,)*,. 

D  /  =  1 


(5.21) 


where  T,t^  is  the  fraction  of  the  total  undamaged  DNA  after  some  exposure 
temperature  and  time.  This  value  is  computed  from  Equation  4.2  with 

/  {T  t  )  = 

Further,  ys  is  the  concentration  of  damaged  total  DNA  which  has  been  repaired  after  time 
tr,  and  yc  is  the  critical  threshold  value  below  which  outgrowth  cannot  occur.  The  domain 
of  integration  D  is  defined  by 

D  =  {(ei,e2,e^):y^(t^,e)  +  Ijj(T,tJ<yJ 


and  fi(ei)  is  the  probability  density  function  for  fitness  due  to  enzyme  that  is 


fM)  = 


r(6r  "h  (3.-1 


nccmjs) 

0, 


-a 


e“  (a-e.y  ,  0<e.<a 

e.  >  a 


(5.22) 
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where  a  =  E^^{T,t^)  .  Now  introducing  the  dimensionless  variables  =  —  and  y.  =  — 

/  a 


from  Section  5.4  into  Equation  5.21  yields 


P,  4,{t,.f)+h(T.t,)<i,  =llJn  /,  r.  — . 

D  i  =  \  ^ 


with 


T{a  +  P)  1 
T{a)T{P)  a 


Yi 


■‘(1-r,) 


p-i 


0, 


0  <  7,  <  1 
7,.  >1 


^  =  {  ( rp  ^2  ’  73  )  :  #5  ’  r )  +  /d  (7^  ’  }  ’ 


^5  represents  the  percentage  of  damaged  DNA  repaired,  represents  the  critical  DNA 
threshold  as  a  percentage  of  the  total  repaired,  and//)f7',tej  is  given  in  Equation  4.2. 

Eigure  5.10  illustrates  the  DNA  repair  for  various  combinations  of  enzymes,  that 
is,  each  curve  represents  the  DNA  repair  within  a  given  spore  from  the  spore  population. 
The  horizontal  axis  represents  the  repair  time  in  seconds.  The  vertical  axis  is  split 
between  the  fraction  of  remaining  DNA  after  an  exposure  time  of  2  seconds  at  a 
temperature  of  SOO^C  ,  foiT,  tj,  and  the  fraction  of  damaged  DNA  repaired,  ^5. 
Approximately  76  %  of  the  DNA  is  still  viable  based  on  the  notional  values  given  in 
Table  4.1,  5.1  and  5.2.  The  red  dashed  line  represents  a  critical  threshold  for  DNA  set  to 
95%  of  total  DNA,  ^c=-95.  If  all  of  the  curves  lie  below  this  critical  threshold  at  an 
arbitrary  chosen  outgrowth  time,  to=150  seconds,  then  the  spore’s  probability  of  survival 
is  zero.  However,  if  the  outgrowth  time  is  chosen  to  be  larger,  say  to=250  seconds,  then 
some  of  the  curves  lie  above  this  critical  threshold  while  others  lie  below.  This  will  lead 
to  a  nonzero  probability  of  survival  for  the  spore. 
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Figure  5.10:  Critical  DNA  Threshold 


The  distribution  of  curves  depicted  in  Figure  5.10  is  dependent  on  the  distribution 
of  enzymes.  As  mentioned  before,  the  combinations  of  viable  enzymes  will  influence  the 
DNA  repair  process.  To  further  explore  this  influence,  consider  different  enzyme 
combinations  that  satisfy 

^5(fr^f)  +  fDiT,te)  =  ^c  (5-23) 

Using  the  nonlinear  model  from  Equation  5.16  with  fixed  exposure  time, 
exposure  temperature,  and  outgrowth  time,  while  each  enzyme  varied  between  zero  and 
one,  it  was  observed  that  different  enzyme  combinations  effected  the  percentage  of 
damaged  DNA  that  is  repaired,  ^5 . 
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Figure  5.11  illustrates  the  different  combinations  of  enzymes  that  will  yield 
particular  threshold  values  as  given  in  Equation  5.23.  The  blue  surface  represents  the 
different  sets  of  enzymes  when  =  .85 .  The  green  surface  represents  the  combination  of 

enzymes  when  =  .90  and  the  red  surface  represents  those  when  =  .95 .  The  domain  of 

integration  for  the  probability  of  kill  would  include  all  combinations  of  enzymes  below 
their  respective  critical  surface.  For  example,  the  domain  of  integration  for  the  probability 
of  kill  when  the  critical  threshold  is  .95  would  include  all  combinations  of  enzymes 
below  the  red  surface  moving  toward  the  origin. 


Figure  5.11:  Enzyme  Combinations  for  Different  Critical  Thresholds 

Based  on  Io=100  nM,  te  =  2  sec,  T  =  350“C,  a=10,  P=3,  to  =  250  sec,  and  values  from 

Table  5.1  and  5.2. 
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After  have  defining  the  probability  of  kill  and  the  domain  of  integration, 
calculations  can  be  made  to  determine  how  outgrowth  time  influences  the  spore’s 
probability  of  survival  after  a  fixed  exposure  temperature  and  time. 

Probability  of  Survival  Calculations. 

Outgrowth  time  is  a  key  factor  for  the  spore’s  probability  of  survival.  For 
example,  after  2  minutes  there  might  not  be  any  visible  vegetative  cells  formed  from  the 
spores.  After  20  minutes  there  might  be  several  visible  vegetative  cells.  This  is  because 
the  spore  has  more  time  to  repair  damaged  DNA.  Therefore,  in  order  to  look  at  the 
spore’s  chance  of  survival,  the  probability  of  survival  will  be  computed  for  three  different 
outgrowth  times  as  well  as  for  two  different  fixed  exposure  temperature  and  time. 

The  probability  of  survival,  Ps,  is  defined  as 

P,=\-P,{^,{f,f)  +  UT,0<0- 

Table  5.4  gives  the  probability  of  survivals  for  different  outgrowth  times.  These 
calculations  were  based  on  4  =  2  seconds,  T  =  300°C,  lo  =  100  nM,  a=  10,  3,  =  -95 

and  the  values  from  Tables  5.1  and  5.2. 

Table  5.4:  Probability  of  Survival  for  Different  Outgrowth  Times  (SOO^C  2  sec) 


Outgrowth  Times 

Probabilitv  of  Survival 

120  seconds 

0% 

240  seconds 

58  % 

300  seconds 

94% 
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From  Table  5.4,  as  the  outgrowth  time  inerease,  the  spore’s  chance  of  survival 
also  increases.  Increased  outgrowth  time  gives  the  spore  additional  time  to  repair 
damaged  DNA  as  expected.  Therefore,  if  the  spore’s  enzymes  are  substantially  damaged, 
the  longer  outgrowth  time  is  crucial  for  the  repair  process  since  the  rate  velocities  will  be 
slow. 

Table  5.5:  Probability  of  Survival  for  Different  Outgrowth  Times  (400'’C  .01  sec) 


Outgrowth  Times 

Probabilitv  of  Survival 

300  seconds 

0% 

540  seconds 

41  % 

720  seconds 

94% 

Table  5.5  gives  the  probability  of  survivals  for  different  outgrowth  times.  These 
calculations  were  based  on  4  =  .01  second,  T  =  dOO^C,  lo  =  100  nM,  a=  10,  =  3,  = 

.95  and  the  values  from  Tables  5.1  and  5.2.  The  increase  in  temperature  allows  pyrolysis 
reactions  to  occur.  Therefore,  the  longer  outgrowth  times  are  required  for  higher 
probability  of  survival. 
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VI.  Conclusions  and  Future  Work 


Conclusions 

In  Chapter  III,  the  thermogravimetric  analysis  data  obtained  from  Dr.  Felker  [18] 
was  examined.  The  purpose  of  this  data  was  to  estimate  the  water  content  within  the 
spore  and  the  temperatures  at  which  this  water  would  become  mobile.  From  these  results 
it  was  concluded  that  the  absorbed  water  becomes  mobile  around  a  peak  of  105°C  and  the 
chemically  bound  water  becomes  available  for  release  at  a  peak  of  290°C.  These  findings 
verify  the  assumptions  made  previously  [19]  who  assumed  that  the  absorbed  water  is 
released  at  a  peak  of  100°C  and  chemically  bound  water  is  released  at  a  peak  of  300°C. 
Additionally,  the  spore’s  weight  was  estimated  to  be  about  40-45%  water.  The  amount  of 
water  within  the  spore  is  important  for  hydrolysis  reactions  to  occur. 

Chapter  IV  examined  the  type  of  damage  occurring  within  the  spore.  Damage  can 
be  caused  by  hydrolysis  and  pyrolysis  reactions.  The  hydrolysis  reactions  can  cause 
damage  to  the  spore’s  DNA  and  enzymes.  This  damage  includes  deamination, 
depurination,  denaturation,  and  depolymerization.  Furthermore,  five  different  DNA  repair 
mechanisms  within  the  spore  were  discussed  with  emphasis  on  BER  as  the  repair 
mechanism  responsible  for  spore  outgrowth.  It  was  assumed  that  a  spore  taken  from  a 
spore  population  begins  with  100%  healthy  DNA  information  and  that  the  DNA 
information  diminishes  exponentially  based  on  rate  coefficients  corresponding  to 
hydrolysis  and  pyrolysis  reactions.  Following  this  further,  the  initial  fitness  distribution 
of  enzymes  developed  by  Captain  Knight  which  followed  a  Gaussian  distribution  was 


78 


enhanced  by  using  a  scaled  beta  distribution  that  corrects  the  concentration  of  enzymes 
within  the  spore  to  fall  between  zero  and  some  maximum  value.  Using  both  the  Gaussian 
and  scaled  beta  distributions,  probability  of  survival  models  were  developed  to  examine 
the  spore’s  chance  of  survival.  These  probability  of  survival  models  were  based  on 
knowledge  of  a  critical  enzyme  threshold,  that  is,  the  value  below  which  the  spore  does 
not  contain  the  sufficient  concentration  of  enzymes  for  the  DNA  repair  process.  Using  the 
notational  values,  it  was  concluded  that  when  heating  the  spore  population  at  a 
temperature  of  225°C  during  dry  heating  and  a  threshold  of  0.10,  the  probability  of 
survival  was  estimated  to  drop  off  faster  when  using  the  beta  representation.  This  is 
because  the  Gaussian  distribution  was  not  bounded  above  as  was  the  beta  distribution. 
However,  because  the  enzyme  critical  threshold  is  unknown  and  the  model  is  compared 
to  actual  data  collected  after  outgrowth,  there  is  a  need  for  an  outgrowth  model. 

The  spore  will  break  dormancy  once  placed  in  a  nutrient  rich  environment.  Once 
this  begins,  the  spore  moves  into  germination.  Germination  occurs  when  the  spore’s 
metabolic  activity  resumes  and  the  spore  begins  repairing  the  damaged  DNA.  It  was 
asserted  that  the  BER  mechanisms  will  be  necessary  to  repair  the  damaged  DNA.  This 
assertion  was  based  on  the  type  of  damage  caused  within  the  spore  from  hydrolysis 
reactions.  Quantitative  models  exist  for  this  DNA  repair  process  within  humans.  Since 
these  repair  process  are  conserved  from  bacteria  to  humans,  these  models  were  mimicked 
to  provide  a  qualitative  understanding  of  the  process  within  the  spore.  Since  the  focus 
was  on  DNA  information  content  instead  of  concentration,  a  dimensionless  variable, 
y-  =  ^/(jwas  introduced  as  was  the  initial  state  of  DNA  within  the  repair  model.  The 
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initial  state  of  DNA  plays  a  major  role  in  the  nonlinear  DNA  BER  model.  It  was 
concluded  that  if  ,^./q  « ,  then  the  DNA  BER  model  could  be  linearized  and  an 

anal3^ical  solution  was  found.  Also,  by  using  the  linear  model,  the  exposure  time 
influence  was  shown  to  appear  only  in  the  exponential  terms.  This  means  that  the 
exposure  time  will  only  change  the  rate  at  which  the  repair  velocity  equilibrates.  When 
^/q  »  ,  then  the  DNA  BER  model  must  remain  nonlinear  and  a  numerical  method 

must  be  used  to  find  the  solution.  Additionally,  since  the  concentration  of  enzyme’s 
initial  distribution  was  based  on  a  beta  distribution,  additional  assumptions  were  made 
about  the  maximum  concentration  for  each  type  of  enzyme  present  during  the  repair 
process.  If  the  spore’s  DNA  is  repaired,  the  spore  then  moves  into  the  last  stage  which  is 
outgrowth.  The  time  at  which  survival  is  assessed  plays  a  crucial  role  in  the 
determination  of  the  spore’s  survival.  Using  this  fact  and  the  enzyme  influence  on  the 
repair  process,  a  probability  of  survival  model  was  developed.  Unlike  before,  this  model 
is  not  based  on  critical  threshold  values  for  the  enzymes.  Instead,  the  probability  of 
survival  model  is  based  on  a  threshold  value  for  the  DNA  information.  Since  the  enzymes 
as  well  as  the  DNA  are  damaged  during  the  thermal  inactivation,  they  influence  the  DNA 
repair  process  which  in  turn  influences  the  probability  of  survival.  Using  the  probability 
of  kill  model  with  the  notational  values,  it  was  concluded  that  with  a  fix  exposure  time  of 
2  seconds,  a  temperature  of  300°C,  and  an  initial  state  of  DNA,  Iq  =  100  nM,  if  given  an 
outgrowth  time  of  120  seconds,  the  spore  has  a  0  %  chance  of  survival.  With  an 
outgrowth  time  of  240  seconds,  the  spore’s  chance  of  survival  increases  to  58  %  and  with 
an  outgrowth  time  of  300  seconds,  the  spore’s  chance  of  survival  increase  to  94  %. 
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Future  Work 


The  Gaussian  and  beta  distributions  were  looked  at  for  the  initial  enzyme 
distribution.  However,  additional  distributions  could  be  considered.  The  DNA  BER 
model  was  based  on  human  DNA  and  even  though,  there  is  comparison  between  bacteria 
and  human  DNA  repair  process,  possible  refinement  in  the  model  should  be  considered 
for  future  work.  Moreover,  to  execute  the  model,  kinetics  data  found  in  literature  based 
on  human  DNA  was  used.  It  would  be  advisable  to  compare  the  model  to  actual 
experimental  data  from  the  spore.  In  addition,  when  developing  the  probability  of 
survival  model  in  Chapter  V  only  the  nonlinear  model  was  considered.  The  linear  model 
might  provide  additional  insight  into  the  spore’s  chance  of  survival.  Further,  assumptions 
about  independence  for  the  enzymes  were  made.  This  assumption  may  not  hold  true, 
therefore  a  joint  distribution  will  need  to  be  consider.  Lastly,  this  research  only 
considered  hydrolysis  reactions  within  the  spore.  Additional  work  should  be  done  that 
looks  at  the  affects  of  the  physical  and  chemical  environment  on  the  spore.  Along  with 
various  chemical  and  biologic  processes  that  can  affect  the  probability  of  survival  of  a 
B.a  spore  after  heat  treatment. 
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